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TARGET  DISCRIMINATION/CLASSIFICATION  RADAR 


1.0  PURPOSE  OF  STUDY 

It  was  proposed  to  investigate  the  use  of  Low  Probability  of  Detection/Super  Wide  Band  (LPI/SWB)  ra¬ 
dar  systems  operating  from  airborne  platforms  for  target  discrimination/classification.  These  LPI/SWB  radar 
systems  would  have  a  transmitted  percent  bandwidth  on  the  order  of  1%  to  10%,  and  they  would  make  use 
of  pulse  compression  techniques.  The  radar  envisioned  would  obtain  its  target  discrimination/classification 
capability  by  exploiting  the  dissimilarities  that  exist  between  the  scattering  characteristics  of  various  radar 
targets  of  interest.  In  other  words,  the  radar  would  have  a  library  of  target  specific  signals  that  it  would  be 
capable  of  transmitting  and  processing.  Each  signal  in  the  library  would  be  constructed  in  such  a  manner 
as  to  excite  a  large  reflection  from  a  specific  target  while  exciting  a  minimal  reflection  from  other  targets. 
Therefore,  when  a  target  is  detected,  the  entire  library  could  be  used  to  interrogate  the  target,  and  the  signal 
that  excited  the  largest  response  would  have  the  highest  probability  of  being  the  target  associated  with  the 
target  specific  signal. 

The  system  would  achieve  its  LPI  characteristics  by  passing  the  target  specific  waveforms  generated  by 
the  target  library  through  a  pulse  expansion  filter.  The  expanded  waveform  would  then  be  up-converted  and 
transmitted.  The  expanded  target  specific  signal  would  then  excite  the  target  and  the  reflected  energy  would 
be  received  and  down-converted.  The  down-converted  signal  would  then  be  passed  through  a  pulse-com¬ 
pression  filter  (which  is  the  time  reversal  of  the  pulse-expansion  filter)  to  obtain  the  minimally  corrupted 
response  of  the  target  to  the  target  specific  signal.  The  system  would  be  similar  to  a  conventional  pulse  com¬ 
pression  system,  and  it  differs  only  in  that  the  standard  pulse  compression  system  excites  the  pulse  expansion 
filter  with  an  impulse  and  the  LPI/SWB  system  excites  the  pulse  expansion  filter  with  the  target  specific  sig¬ 
nal. 

It  was  anticipated  that  this  investigation  would  be  achieved  in  three  steps: 

1)  Various  LPI/SWB  radar  implementations  would  be  rigorously  investigated  to  determine  which 
method  has  the  most  promise  in  the  application  described  above. 

2)  The  scattering  characteristics  of  various  target  geometries  of  interest  would  be  investigated  to 
determine  exactly  how  they  may  be  exploited. 

3)  The  information  learned  in  the  studies  outlined  above  would  be  combined  to  determine  guidelines 
that  can  be  used  in  the  design  of  specific  radar  systems  to  accomplish  the  target  discrimination/ 
classification  task  described  above. 

Task  1  involved  an  investigation  into  various  types  of  pulse  compression  radar  systems  (PSK,  FSK  and 
Hybrid  FSK/PSK)  to  determine  their  limitations  and  beneficial  qualities.  Some  of  these  qualities  investi¬ 
gated  were  achievable  bandwidth,  low  probability  of  intercept  capability,  and  the  spectral  purity  of  the  pulse 
expansion  filter.  Task  2  involved  the  development  of  accurate  computer  models  of  generic  targets  for  genera¬ 
tion  of  target  scattering  characteristics.  This  involved  the  use  of  existing  exact  solution  techniques  that  could 
be  used  for  resonant  region  analysis  (or  the  region  at  which  the  target  size  is  on  the  order  of  a  wavelength), 
and  the  development  of  approximate  solution  programs  for  optical  region  analysis  (or  the  frequency  region 
where  the  target  extends  for  many  wavelengths).  Task  3  involved  using  the  tools  developed  in  Tasks  1  and 
2  to  determine  the  optimum  implementation  of  a  LPI/SWB  radar.  This  was  basically  an  analysis  of  model 
systems  to  determine  performance  levels  as  target  discrimination/classification  systems. 
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2.0  SUMMARY  AND  CONCLUSIONS  OF  RESEARCH 


The  purpose  of  this  study  was  to  develop  a  signal  processing  scheme  that  would  combine  the  separate 
research  efforts  of  low  probability  of  intercept  signaling,  low  altitude  low  observable  target  tracking,  and  non- 
cooperative  target  recognition  into  a  single  research  effort.  The  end  result  of  this  fusion  is  a  new  radar  design 
concept  designated  adaptive  FSK/PSK  signaling.  An  adaptive  FSK/PSK  radar  system  makes  use  of  matched 
FSK/PSK  signals  in  conjunction  with  a  gradient  descent  adaptive  signal  design  algorithm.  Matched  FSK/ 
PSK  signals  are  ultra  wideband  pseudorandom  multiple  frequency  bi-phase  modulated  signal  sequences  that 
have  spectral  content  specifically  matched  to  the  target  under  track.  This  use  of  target-specific  information 
in  the  design  of  radar  signals  is  the  salient  point  of  this  research  effort.  It  allows  for  more  efficient  use  of  signal 
power  and  results  in  a  performance  improvement  over  existing  radar  technologies  that  are  matched  to  the 
transmit  signal  only.  Since  the  spectral  content  of  the  target  is  not  known  a  priori,  the  target  range  profile 
must  be  learned  by  the  radar.  This  matching  of  the  signal  to  the  target  is  accomplished  through  the  use  of 
an  iterative  (adaptive)  algorithm  that  is  based  on  gradient  descent.  Once  the  radar  has  adapted  to  the  target 
(learned  the  target),  the  range  profile  of  the  target  is  known  and  target  identification  can  take  place  based  on 
the  learned  range  profile.  Therefore,  an  adaptive  FSK/PSK  radar  system  can  be  thought  of  as  a  learning  sys¬ 
tem  as  opposed  to  a  measmement  system.  In  other  words,  the  radar  measurement  is  improved  as  a  result  of 
the  learning  process. 

The  pseudorandom  sequence  structure  of  the  matched  FSK/PSK  signal,  along  with  its  high  time-band¬ 
width  product,  makes  the  radar  inherently  difficult  to  detect.  Therefore,  the  signals  can  be  classified  as  low 
probability  of  intercept.  Additionally,  since  the  transmitter  power  is  used  in  a  most  effective  manner  (signal 
efficiency  is  maximized  by  target  matching),  lower  peak  power  signals  can  be  utilized  by  the  radar.  This  re¬ 
sults  in  an  additional  improvement  in  LPI  performance. 

In  an  additive  noise  environment,  matched  FSK/PSK  signaling  provides  a  signal  processing  gain  over 
a  baseline  radar  that  utilizes  a  signal  with  an  impulse-like  ambiguity  function.  This  performance  improve¬ 
ment  results  from  the  fact  that  the  signal  power  is  utilized  in  an  efficient  manner  and  the  target  reflection  is 
enhanced  without  enhancing  the  additive  noise.  The  level  of  signal  processing  gain  is  dependent  upon  the 
spectral  complexity  of  the  target.  If  the  target  reflects  all  frequencies  equally  well,  little  processing  gain  is 
seen.  On  the  other  hand,  if  the  target  reflectivity  is  concentrated  in  a  few  isolated  frequency  bands,  then  the 
processing  gain  can  be  large  (on  the  order  of  the  square  of  the  number  of  frequencies  used  in  the  signal). 

In  a  ground  clutter  limited  environment,  matched  FSK/PSK  signaling  performs  better  than  conventional 
radar  technologies  that  do  not  take  into  account  any  target  specific  information.  This  occurs  because  the 
matched  FSK/PSK  radar  system  is  capable  of  achieving  a  theoretical  target  enhancement  of  between  1  and 
K^  while  the  integrated  sidelobe  level  is  increased  by  a  factor  of  between  1  and  K  where  K  is  the  number  of 
frequencies  used  in  the  construction  of  the  signal.  Therefore,  in  a  clutter  limited  environment,  adaptive  FSK/ 
PSK  radar  systems  can  achieve  performance  improvement  factors  of  between  1  and  K.  It  should  be  noted 
that  a  matched  FSK/PSK  radar  system  and  a  conventional  radar  system  using  a  matched  filter  post-processor 
are  able  to  achieve  identical  performance  in  a  clutter  limited  environment. 

Analysis  of  predicted  real  target  signature  data  shows  that  the  use  of  target-specific  information  in  the 
design  of  the  radar  signal  typically  yields  average  signal  to  noise  ratio  improvements  of  between  3  dB  and 
6  dB  depending  on  the  target.  In  a  clutter  limited  environment,  typieal  improvements  are  between  2  dB  and 
5  dB .  These  improvements  are  closer  to  the  low  end  of  the  theoretically  possible  performance  improvements. 
This  is  result  of  the  uniform  nature  of  optical  region  target  spectra.  The  gains,  however,  are  significant  and 
provide  a  significant  performance  increase.  It  should  be  pointed  out  that  in  the  resonance  region,  greater  gains 
are  expected  since  resonance  region  spectra  are  typically  impulse-like. 

The  adaptation  algorithm  derived  in  the  research  (See  Chapter  VIII  of  the  Appendix)  provides  a  method 
of  not  only  designing  the  optimum  signal  to  excite  the  maximum  response  from  the  target,  but  also  for  mea¬ 
suring  the  best  mean  squared  estimate  of  the  target  range  profile.  Through  computer  simulations,  it  was 
shown  that  the  adaptation  algorithm  removes  additive  noise.  The  level  of  noise  removal  is  dependent  on  the 
value  of  the  learning  parameter.  The  smaller  the  learning  parameter  (slower  the  learning  rate),  the  better  the 
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mean  squared  estimate.  The  simulations  also  show  that  the  random  fluctuations  of  the  target  range  profile 
due  to  random  range  motion  and  random  incidence  angle  motion  increase  the  mean  squared  error  of  the  esti¬ 
mate.  If  the  learning  rate  is  high  enough  for  the  adaptive  algorithm  to  learn  the  target  range  profile  before 
the  profile  decorrelates,  then  the  errors  are  small.  If  on  the  other  hand,  the  target  decorrelates  more  rapidly 
than  the  algorithm  can  learn,  range  profile  estimation  errors  are  increased.  This  results  in  decreased  target 
enhancement  due  to  matched  FSK/PSK  signaling.  It  should  be  pointed  out  that  the  performance  metric  used 
in  this  study  (the  mean  squared  range  profile  error)  is  biased  towards  stationary  (nonfluctuating)  targets.  It 
is  difficult  to  ascertain  the  deleterious  effects  of  range  profile  fluctuations  on  an  actual  target  identification 
algorithm.  The  effects  may  be  insignificant  since  the  adaptation  algorithm  determines  a  best  mean  squared 
fit  to  the  second  order  statistics  of  the  random  target  fluctuation  process  within  the  learning  period  of  the  algo¬ 
rithm.  It  is  of  interest  to  note  that  a  linear  averaging  scheme  would  yield  a  zero  result  for  range  profile  estima¬ 
tion  since  the  phase  angle  of  the  returns  is  random  from  measurement  to  measurement. 

In  conclusion,  adaptive  FSK/PSK  signaling  radar  is  an  excellent  platform  for  fulfilling  the  future  mission 
of  radar:  low  probability  of  intercept  noncooperative  target  recognition.  It  utilizes  low  probability  of  inter¬ 
cept  signals  in  a  spectrally  adaptive  manner.  The  ability  of  the  radar  to  adapt  its  spectrum  makes  it  possible 
to  enhance  targets  thus  improving  system  performance  in  additive  noise  environments  and  ground  clutter  lim¬ 
ited  environments.  Finally,  the  adaptation  scheme  derived  in  this  study  is  capable  of  extracting  the  range 
profile  of  the  target  under  track.  This  is  important  for  target  identification  since  it  is  the  range  profile  that 
is  the  feature  vector  used  by  most  target  identification  algorithms. 

In  the  work  performed,  a  series  of  equations  were  derived  to  model  the  Radar  Performance  for  Target 
Return,  Signal-to-Clutter  and  Bandwidth  Loss  for  three  cases.  These  three  cases  are  a  Baseline  Radar  (Cor¬ 
relation)  Receiver,  a  Matched  Filter  Radar  Receiver,  and  a  Matched  FSK/PSK  Radar  Receiver.  The  Table 
below  summarizes  the  equations  with  equation  numbers  referencing  Chapter  VII  of  the  Appendix. 

Performance  estimates  for  specific  target  scenarios  are  contained  in  Chapter  VII  for  Appendbc,  Table  7.1, 
page  91. 

Summary  of  Radar  Performance  Equations. 


Baseline 

Matched  Filter 

Matched  FSK/PSK 

Target  Return 
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Equation  (7.18) 

Signal  to  Clutter 

Eq 
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2 
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Equation  (7.22) 

Bandwidth  Loss 
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Equation  (S.14) 
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APPENDIX 


The  appendix  is  a  reproduction  ofthe  dissertation  that  resulted  from  the  research  for  this  project.  Publica¬ 
tions  that  occurred  during  the  course  of  the  research  are  listed  in  Section  3.0.  The  appendix  is  a  self  contained 
document  that  gives  full  background  details  and  results.  It  is  the  desire  of  the  authors  that  the  work  be  readable 
and  useable  to  those  who  are  interested  in  this  interesting  and  vital  topic. 
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CHAPTER  I 


INTRODUCTION 

Radar  was  first  used  as  a  military  sensor  system  during  World  War  n  when  it 
came  into  prominence  during  the  air  war  over  Britain.  These  early  radar  systems  were 
able  to  locate  the  direction  and  distance  (azimuth  and  range)  of  attacking  waves  of 
German  bombers  while  they  were  stiU  weU  beyond  visual  range.  This  location 
information  made  it  possible  to  deploy  the  outnumbered  Royal  Air  Force  fighter  wings 
in  a  more  effective  manner  and  thus,  the  Battle  of  Britain  was  won.  By  the  end  of  the 
war,  the  state-of-the-art  in  radar  had  advanced  significantly  and  included  fire  control 
radar  for  anti-aircraft  artillery.  These  radar  systems  were  able  to  automatically  track 
the  range,  range-rate,  azimuth,  and  elevation  of  a  target. 

During  the  late  1940’s  and  1950’s,  the  role  of  radar  technology  remained 
basically  unchanged.  Developments  in  radar  were  primarily  devoted  to  codifying  radar 
design  principles  [1]  and  expanding  the  underlying  mathematical  theory  of  radar  [2]. 
During  the  1960’s  and  1970’s,  the  new  technologies  of  solid  state  devices  and  digital 
computing  were  implemented  into  radar  systems  making  them  smaller,  more  power 
efficient,  and  basically  more  capable.  StiU,  technology  had  not  advanced  to  the  extent 
that  the  role  of  World  War  n  era  radar  would  be  changed  significantly  change. 
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The  1970’s  and  1980’s  saw  the  advent  of  two  new  radar  design  requirements: 
tracking  of  stealth  (or  low  observable)  aircraft  and  anti-radiation  missiles.  Stealth 
technology  drove  research  into  the  area  of  radar  signal  designs  capable  of  tracking  low 
radar  cross  section  targets  buried  in  ground  clutter  returns.  The  goal  of  these  signal 
designs  was  waveforms  capable  of  rejecting  the  range  extended  ground  clutter  return 
[3,  4,  5,  6]  primarily  by  lowering  the  range  sidelobes  of  the  radar  signal  ambiguity 
function.  These  clutter  suppression  efforts  assumed  that  the  radar  target  of  interest 
was  localized  to  a  point-like  range  extent  (point  target  assumption)  and  therefore, 
reflected  all  waveforms  equally  well.  Anti-radiation  missiles  drove  research  into  the 
area  of  signal  designs  that  made  it  more  difficult  for  a  radar  warning  receiver  to  detect 
the  existence  of  the  radar  signal  (low  probability  of  intercept  signaling)  [7].  In  other 
words,  the  radar  signals  themselves  were  forced  to  become  stealthy.  This  led  to  pulse 
compression  signals  that  utilized  random  frequency  and  phase  coding  to  provide  a 
signal  processing  gain  to  the  radar.  This  signal  processing  gain  was  typically 
measured  by  the  time-bandwidth  product  of  the  waveform  [7,  8,  9]. 

The  late  1980’s  and  1990’s  have  seen  an  interest  in  expanding  the  role  of 
radar.  Instead  of  radar  simply  being  a  detection  device,  radar  systems  capable  of 
identifying  targets  [10]  are  of  interest.  In  other  words,  modern  radar  systems  will  be 
required  to  extract  more  information  from  targets.  Therefore,  there  is  an  increased 
interest  in  so-called  ultra  wideband  radar  concepts  [10].  Ultra  wideband  radar 
concepts  can  be  classified  into  two  types  of  radar  signaling  regimes  [11]:  impulse 
radar  for  resonant  frequency  extraction  and  high  range  resolution  waveform  signaling. 
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Resonant  region  radars  try  to  identify  targets  based  on  the  extraction  of  the 
natural  resonant  frequencies  excited  by  a  broadband  pulse  of  electromagnetic  radiation 
[12,  13,  14,  15,  16].  These  resonant  frequencies  are  target-specific  and  are  therefore 
good  classification  features.  The  problem  with  this  technique  is  that  the  radars  are 
required  to  operate  in  the  VHP  and  low  UHF  frequency  bands  and  therefore,  the  radar 
systems  are  very  large  and  unwieldy.  Additionally,  the  percent  bandwidth  of  the 
waveforms  used  must  be  very  high  (multi-octave  signaling  desirable)  and  Doppler 
processing  becomes  ineffective  for  ground  clutter  removal. 

High  range  resolution  radar  systems  identify  radar  targets  based  on  the  location 
of  scattering  centers  [17,  18].  In  other  words,  the  range  resolution  of  the  radar  is  such 
that  the  individual  scattering  components  of  the  target  can  be  resolved.  These  radar 
systems  operate  in  typical  radar  frequency  bands  (above  1  GHz),  but  use  higher  than 
typical  signal  bandwidths.  High  range  resolution  radar  systems  have  the  advantage  of 
small  size,  narrow  antenna  beamwidths,  and  a  higher  degree  of  Doppler  precision. 

The  disadvantage  of  high  range  resolution  radar  systems  is  that  the  range  profiles  are 
highly  dependent  upon  target  aspect  angle  and  a  large  amount  of  data  is  required  to 
describe  each  individual  target  in  a  target  database  [17,  18]. 

1.1.  Combining  Goals 

Future  radar  systems  will  use  low  probability  of  intercept  (LPI)  waveforms  to 
track  and  identify  low  observable  targets.  Therefore,  the  development  of  radar  target 
identification  techniques  must  be  accomplished  within  the  confines  of  the  requirements 
on  low  observable  target  tracking  and  signal  detectability  by  hostile  interceptors.  In 
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the  past,  research  in  these  three  areas  has  progressed  in  separate  directions.  Research 
into  low  observable  target  tracking  has  been  based  on  the  assumption  that  radar  targets 
are  point  targets  (single  point  scatterer)  while  ground  clutter  is  range  extended. 
Therefore,  design  techniques  have  focused  on  designing  waveforms  that  are  matched 
to  point  targets  yet  reject  distributed  clutter  targets.  Target  identification  algorithms, 
on  the  other  hand,  make  use  of  the  range  extended  nature  of  radar  targets  [17,  18] 

(i.e.,  each  target  scatters  radiation  differently).  A  question  therefore  arises;  should  low 
observable  target  tracking  research  focus  on  enhancing  the  return  from  distributed 
radar  targets  by  making  use  of  target-specific  information  as  well  as  suppressing 
clutter?  In  short,  is  it  possible  for  the  range  extended  target-specific  information  to  be 
utilized  to  the  radar’s  advantage? 

t 

1.2.  Scope  of  Study 

The  purpose  of  this  study  is  to  combine  the  three  different  research  efforts 
related  to  low  observable  target  tracking  in  clutter,  LPI  signal  design,  and  target 
identification  into  a  single  research  effort.  In  this  study,  a  pulse  compression 
technique  designated  "matched  FSK/PSK"  [21]  is  introduced.  Unlike  conventional 
signaling  techniques,  matched  FSK/PSK  signaling  makes  it  possible  to  not  only 
suppress  clutter  returns,  but  also  to  enhance  the  target  return.  The  target  return 
enhancement  is  accomplished  by  making  use  of  target-specific  information  in  the 
signal  design  process.  Since  target  enhancement  requires  a  priori  knowledge  about  the 
range  profile  of  the  particular  target,  an  iterative  signal  design  technique  [22]  is 
introduced.  This  is  a  nonlinear  design  technique  that  effectively  extracts  the  target’s 
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range  profile  over  time  and  can  be  thought  of  as  a  measurement  technique.  The 
extracted  target  range  profile  is  a  usable  feature  for  target  identification  purposes. 

When  matched  FSK/PSK  is  used  in  conjunction  with  the  iterative  design  algorithm, 
the  resulting  system  is  designated  "adaptive  FSK/PSK".  This  study  is  theoretical  in 
scope  since  a  treatment  of  hardware  implementation  issues  would  tend  to  mask  the 
fundamental  principles  involved.  Realistic  target  modeling  is  used,  however,  to 
acquaint  the  reader  with  typical  achievable  performance  improvements. 

1.3.  Organization  of  Study 

This  study  is  organized  in  the  following  manner.  Some  fundamentals  of  radar 
and  radar  signal  processing  are  presented  briefly  in  Chapter  n.  The  concepts  of  pulse 
compression,  radar-interceptor  advantage,  range  smearing,  ground  clutter  analysis  with 
the  cross-ambiguity  function,  as  well  as  radar  scattering  regimes  are  discussed  in 
general  terms  in  order  to  standardize  terminology  and  to  develop  a  basis  for  presenting 
and  solving  the  problem.  In  Chapter  m,  the  cross-ambiguity  function  [20]  is 
discussed  in  general  with  respect  to  radar-target  interaction,  noise  removal  and 
resolution,  and  specifically  with  respect  to  FSK/PSK  signaling.  The  cross-ambiguity 
function  of  the  FSK/PSK  waveform  is  analyzed  probabilistically  in  Chapter  IV.  This 
analysis  includes  a  mean  and  mean  squared  analysis  of  the  mainlobe  of  the  cross¬ 
ambiguity  function,  mean  and  mean  squared  analysis  of  the  cross-ambiguity  function 
sidelobes,  and  evaluation  of  the  integrated  mean  square  sidelobe  level  of  the  cross¬ 
ambiguity  function  (a  baseline  measure  of  clutter  rejection).  The  power  spectral 
density  and  equivalent  time-bandwidth  product  of  an  FSK/PSK  radar  signal  are 
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derived  in  Chapter  V.  In  Chapter  VI,  the  probabilistic  analysis  of  Chapter  IV  is  used 
to  develop  a  design  methodology  for  matched  FSK/PSK  radar  based  on  a  known 
arbitrary  range  profile.  The  performance  gain  of  matched  FSK/PSK  over  a 
conventional  signaling  radar  system  that  utilizes  no  target-specific  information  is 
determined  in  Chapter  VH.  The  performance  metrics  used  in  the  study  are  the  target 
enhancement,  signal  to  noise  ratio,  and  signal  to  clutter  ratio  of  matched  FSK/PSK 
relative  to  a  baseline  conventional  radar.  These  performance  gains  are  calculated  for 
both  theoretical  target  scattering  distributions  and  for  some  predicted  target  data.  An 
iterative  signal  design  technique  is  derived  in  Chapter  VUI.  This  technique  is  derived 
based  on  gradient  descent  optimization  of  the  mean  squared  reflected  voltage.  This 
gradient  descent  adaptive  signal  design  scheme  is  also  studied  parametrically  in 
Chapter  VIH.  This  parametric  study  illustrates  the  effect  of  additive  white  noise, 
target  range  motion,  and  target  angular  motion.  Concluding  remarks  and  suggestions 
for  future  work  are  discussed  in  Chapter  DC. 

1  A.  References 

1.  Radiation  Laboratory  Reports 

2.  P.  M.  Woodward,  Probability  and  Information  Theory,  with  Applications  to 
Radar.  Pergamon  Press,  New  York,  1953. 

3.  A.  W.  Rihaczek  and  R.  L.  Mitohell,  "Radar  Waveforms  for  Suppression  of 
Extended  Clutter",  IEEE  Transactions  on  Aerospace  and  Electronic  Systems. 
May  1967,  pp.  656-663. 

4.  D.  F.  DeLong,  Jr.  and  E.  M.  Hostetter,  "On  the  Design  of  Optimum  Radar 
Waveforms  for  Clutter  Rejection",  IEEE  Transactions  on  Aerospace  and 
Electronic  Systems.  January  1967,  pp.  138-141. 


7 


5.  M.  Ares,  "Optimum  Burst  Waveforms  for  Detection  of  Targets  in  Uniform 
Range-Extended  Clutter",  IEEE  Transactions  on  Aerospace  and  Electronics 
Systems,  January  1967,  pp.  138-141. 

6.  T.  K.  Bhattacharya  and  P.  R.  Mahapatra,  "A  Versatile  and  Powerful  Method  of 
Designing  Optimum  Signals  Against  Arbitrary  Range-Doppler  Clutter",  1989 
International  Radar  Conference. 

7.  G.  Schrick  and  R.  G.  Wiley,  "Interception  of  LPI  Radar  Signals",  1990  IEEE 
International  Radar  Conference,  pp.  108-111. 

8.  D.  C.  Schleher,  "Low  Probability  of  Intercept  Radar",  1985  IEEE  International 
Radar  Conference,  pp.  346-349. 

9.  S.  L.  Johnston,  "Radar  Operability/Survivability  in  the  1990’s",  1989 
International  Radar  Conference. 

10.  M.  Skolnik,  "Status  of  Ultra  Wideband  (UWB)  Radar  and  its  Technology", 
Proceedings  of  the  1992  IEEE  Antennas  and  Propagation  Society  International 
Symposium.  Vol.  3.  pp.  1224-1227. 

11.  L.  S.  Riggs  and  C.  R.  Smith,  "Bayesian  Probability  Theory  Applied  to  the 
Problem  Radar  Target  Discrimination",  Proceedings  of  the  1992  IEEE  Antennas 
and  Propagation  Society  International  Symposium.  Vol.  3.  pp.  1246-1249. 

12.  A.  J.  Bemi,  "Target  Identification  by  Natural  Resonance  Estimation",  IEEE 
Transactions  on  Aerospace  and  Electronics  Systems.  March  1975,  pp.  147-154. 

13.  A.  A.  Ksienski,  Y.T.  Lin,  and  L.  J.  White,  "Low-Frequency  Approach  to  Target 
Identification",  Proceedings  of  the  IEEE.  December  1975,  pp.  1651-1660. 

14.  C.  W.  Chuang  and  D.  L.  Moffatt,  "Natural  Resonances  of  Radar  Targets  Via 
Prony’s  Method  and  Target  Discrimination",  IEEE  Transactions  on  Aerospace 
and  Electronics  Systems.  September  1976,  pp.  583-589. 

15.  E.  W.  Kennaugh,  "The  K-Pulse  Concept",  IEEE  Transactions  on  Antennas  and 
Propagation.  March  1981,  pp.  327-331. 

16.  E.  K.  Rothwell,  D.  P.  Nyquist,  K.  M.  Chen,  and  B.  Drachman,  "Radar  Target 
Discrimination  Using  the  Extinction-Pulse  Technique",  IEEE  Transactions  on 
Antennas  and  Propagation.  September  1985,  pp.  929-936. 


8 


17.  S.  Hudson  and  D.  Psaltis,  "Correlation  Filters  for  Aircraft  Identification  From 
Radar  Range  Profiles",  IEEE  Transactions  on  Aerospace  Electronics  Systems, 
July  1993,  pp.  741-748. 

18.  H.  J.  Li  and  S.  H.  Yang,  "Using  Range  Profiles  as  Feature  Vectors  to  Identify 
Aerospace  Objects",  IEEE  Transactions  on  Antennas  and  Propagation.  March 
1993,  pp.  261-268. 

19.  E.  J.  Kelly  and  R.  P.  Wishner,  "Matched-Filter  Theory  for  High-Velocity, 
Accelerating  Targets",  IEEE  Transactions  on  Military  Electronics.  January 
1965,  pp.  56-69. 

20.  A.  W.  Rihaczek,  Principles  of  High  Resolution  Radar.  Peninsula  Publishing, 
Los  Altos,  CA,  1985. 

21.  B.  J.  Skinner,  J.  P.  Donohoe,  and  F.  M.  Ingels,  "Matched  FSK/PSK  Radar", 
Proceedings  of  the  1994  National  Radar  Conference,  pp.  251-255. 


22. 


B.  J.  Skinner,  J.  P.  Donohoe,  and  F.  M.  Ingels,  "Gradient  Descent  Method  for 
Designing  Optimum  Radar  Signals",  accepted  for  publication  Proceedings  of 
the  1995  International  Radar  Conference. 


CHAPTER  n 


RADAR  FUNDAMENTALS 

In  this  chapter,  several  fundamentals  of  radar  systems  analysis  are  discussed 
primarily  to  introduce  notation  and  for  proper  problem  formulation.  These 
fundamental  topics  are:  1)  signal  to  noise  ratio  (SNR)  at  the  radar  receiver  and  at  the 
intercept  receiver  (receiver-interceptor  advantage),  2)  pulse  compression,  3)  range 
resolution,  4)  velocity  or  Doppler  resolution,  5)  range  smearing,  6)  range-doppler 
trade-offs,  and  7)  ground  clutter  analysis  with  the  cross-ambiguity  function. 
Additionally,  an  overview  of  radar  scattering  fundamentals,  as  related  to  target 
identification,  is  given.  These  limits  and  bounds  are  not  derived  in  detail,  only 
discussed  in  relation  to  the  problem  of  LPI  tracking  and  identification  of  low 
observable  radar  targets. 

Figure  2.1  graphically  shows  all  components  of  the  problem  that  are  points  of 
interest  in  this  study.  A  waveform  y(t)  is  transmitted  ft'om  the  radar  anterma.  A 
portion  of  the  energy  is  reflected  by  the  target  while  another  portion  is  reflected  by  the 
surface  of  the  earth.  The  energy  reflected  ft'om  the  earth  is  referred  to  as  clutter  as  it 
is  an  unwanted  return.  The  received  signal  and  clutter  is  also  contaminated  with 
additive  noise  from  various  sources  such  as  thermal  noise,  the  sun,  other  radars  and 
communication  equipment  in  the  area,  and  hostile  jamming  equipment.  Once  the 
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Figure  2.1.  General  Radar  Block  Diagram  Showing  the  Target  and  a  Clutter 
Contaminant. 

signal  is  received  by  the  antenna,  it  is  passed  to  a  signal  processor  that  processes  the 
waveform  and  separates  out  the  returns  by  range  and  Doppler  (velocity).  Typically  the 
clutter  is  stationary  whereas  the  target  is  moving  toward  or  away  from  the  radar. 
Therefore,  the  target  and  clutter  theoretically  can  be  separated  based  on  relative 
velocity.  In  reality,  however,  measurement  of  range  and  velocity  are  not  perfect,  and 
clutter  can  obscure  the  moving  target.  This  will  be  discussed  in  more  detail  later  in 
this  chapter. 


2.1.  The  Radar  Equation 

The  radar  equation,  which  describes  the  signal  to  noise  ratio  of  the  processed 


received  waveform  is  given  by 
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SNR^ 


P,  G '  o 
(471)^ 


(2.1) 


where  P,  is  the  transmitter  peak  power,  is  the  antenna  gain.  A,  is  the  radar 
wavelength,  a  is  the  effective  radar  cross  section  of  the  target  within  the  radar 
bandwidth,  is  the  noise  power  spectral  density.  Bp  is  the  processor  bandwidth,  and 
R  is  the  range  to  the  target.  A  substitution  of 
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where  is  the  antenna  effective  aperture,  into  (2.1)  yields  the  following  equation  in 
terms  of  the  antenna  aperture  and  wavelength 


SNR  - _ ° _  .  (2.3) 

(471)  A"  A  B 

Therefore,  for  a  fixed  antenna  aperture  A^,  the  received  power  increases  with 
decreasing  wavelength. 

The  signal  to  noise  ratio  for  a  intercept  receiver  on  board  the  target  being 
illuminated  by  the  radar  is  given  by 


SNR. 


P,  G^  G,.  A^ 
5  (  471  )  R^ 


(2.4) 


where  G,  is  the  antenna  gain  of  the  intercept  receiver  and  B^  is  the  radiated  bandwidth 
of  the  radar  signal.  In  Equation  (2.4),  the  noise  levels  are  assumed  to  be  the  same  in 
both  the  intercept  receiver  and  the  radar  receiver  and  the  intercept  receiver  is  assumed 
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to  be  in  the  mainlobe  of  the  radar  antenna.  It  is  also  assumed  that  the  radiated 
bandwidth  of  the  radar  signal  is  much  greater  than  the  bandwidth  of  the  intercept 
receiver  signal  processor.  An  important  quantity  is  the  ratio  of  the  radar  SNR  to  the 
intercept  receiver  SNR.  This  is  known  as  the  radar-interceptor  advantage  a  and  is 
given  by 


a  = 


SNR^ 

Wr. 


f  \ 
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The  radar-interceptor  advantage  is  important  because  it  relates  the  performance  of  the 
radar  to  the  performance  of  the  radar’s  enemy  (the  intercept  receiver).  Signal  to  noise 
ratios  are  used  because  they  typically  describe  the  "bottom  line"  of  system 
performance.  Therefore,  if  the  radar  requires  a  minimum  signal  to  noise  ratio  of 
SNR^i„  to  meet  its  specified  performance  level,  and  the  intercept  receiver  requires  a 
minimum  signal  to  noise  ratio  of  SNRi^„  to  meet  its  specified  performance  criterion, 
the  radar  system  will  have  an  advantage  over  the  intercept  receiver  if 


a  > 


SNR 


SNR,_ 


(2.6) 


(assuming  the  minimum  radar  SNR  has  been  achieved).  Otherwise,  the  intercept 
receiver  is  said  to  have  the  advantage.  Solving  Equation  (2.5)  for  the  range  R^  at 
which  the  radar  achieves  an  advantage  over  the  intercept  receiver  yields 
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where  is  the  minimum  radar-interceptor  advantage  satisfying  Equation  (2.6). 
Therefore,  if  the  minimum  SNR  for  the  radar  can  be  lowered,  the  operable  range  of 
the  radar  can  be  increased.  Notice  that  there  is  a  5  dB/decade  improvement  in  range 
versus  radar  signal  to  noise  ratio  improvement.  Notice  also  that  there  are  three  factors 
available  to  the  radar  designer  in  Equation  (2.7)  that  improve  radar  performance: 
antenna  gain  relative  to  the  intercept  receiver  anteima  gain  G„  radar  radiated 
bandwidth  relative  to  the  radar  receiver  processor  bandwidth  Bp,  and  possibly  the 
effective  radar  cross  section  a  of  the  target.  The  antenna  gain  parameter  is  a  hardware 
parameter  that  is  fixed  by  the  operating  frequency  of  the  radar  along  with  the  physical 
dimensions  of  the  antenna.  The  ratio  of  radiated  bandwidth  to  processor  bandwidth  is 
a  signal  processor  dependent  parameter  that  can  be  improved  by  utilizing  advanced 
signal  processing  schemes  such  as  pulse  compression.  The  effective  radar  cross 
section  of  the  target  is  typically  thought  of  as  a  parameter  that  is  fixed  by  the  target. 

In  this  study,  however,  signal  processing  schemes  will  be  introduced  that  make  use  of 
target-specific  information  in  order  to  increase  the  effective  radar  cross  section  of  the 
target  thus  increasing  the  LPI  performance  of  the  radar.  Therefore,  effective  radar 
cross  section  can  be  thought  of  as  a  design  parameter  available  to  the  radar  designer. 
In  this  study,  only  signal  processing  gains  are  investigated. 


2.2.  Pulse  Compression 


Pulse  compression  is  a  signaling  teclmique  in  which  a  signal  of  duration  and 


bandwidth  B^  is  transmitted  and  the  reflected  waveform  is  compressed  through  a 
matched  filter  into  a  pulse  of  duration  HB^.  This  makes  it  possible  for  a  radar  to 
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transmit  a  long  duration  low  peak  power  signal  that  can  be  compressed  into  a  short 
duration  high  peak  power  signal.  This  is  shown  graphically  in  Figure  2.2.  Another 
interpretation  and  implementation  of  pulse  compression  involves  the  use  of  a 
correlator.  When  a  correlator  is  used  for  despreading,  the  signal  of  bandwidth  is 
coUapsed  into  a  signal  of  bandwidth  l/T^.  This  is  shown  graphically  in  Figure  2.3. 

For  a  pulse  compression  radar,  Equation  (2.5)  becomes 


a  = 


^  1  'l  ('  1  'l 


471 

V  y 


(2.8) 


Since  is  always  greater  than  unity  if  there  is  any  pulse  compression  effect,  the 
radar-interceptor  advantage  is  always  increased.  Pulse  compression  is  desirable 
because  it  allows  a  radar  to  have  the  range  resolution  of  a  short  pulsewidth  radar 
signal  with  the  peak  power  requirements  of  a  long  pulsewidth  radar  signal. 


2.3.  Range  Resolution 

The  range  resolution  of  a  radar  is  defined  to  be  the  minimum  distance  between 
two  point  targets  at  which  the  radar  can  determine  that  there  are  two  targets  and  not 
one.  This  is  shown  graphically  in  Figure  2.4.  The  range  resolution  A/?  of  a  radar 
transmitting  a  signal  with  a  bandwidth  B^  is  given  by  [1,  2] 


AR  = 


c 


2B 


(2.9) 


where  c  is  the  speed  of  light.  Therefore,  as  the  bandwidth  of  the  signal  increases,  the 
range  resolution  increases. 


Figure  2.3.  Correlator  Implementation  of  Pulse  Compression. 
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2.4.  Doppler  Resolution 

Doppler  resolution  (or  precision)  is  the  minimum  difference  in  the  velocity  of 
two  targets  at  which  the  radar  can  determine  that  there  are  two  targets  and  not  one. 
This  is  shown  graphically  in  Figure  2.5.  The  Doppler  resolution  Av  of  a  radar 
transmitting  a  signal  of  duration  is  given  by 

Av  = _ —  (2.10) 

2/o7’, 

Therefore,  the  ability  of  a  target  to  discriminate  (or  separate)  two  targets  moving  at 
two  different  velocities  is  related  directly  to  the  duration  of  the  signal  and  the  carrier 
frequency  of  the  radar. 

2.5.  Range  Smearing 

The  previous  two  identities  would  lead  one  to  believe  that  the  key  to  radar 
performance  would  be  to  have  an  infinite  bandwidth  infinite  duration  signal. 
Physically  this  cannot  occur  because  of  range  smearing.  Range  smearing  occurs  with 
pulse  compression  radars  when  a  target  moves  an  appreciable  portion  of  a  range 
resolution  cell  during  the  signal  dmation.  To  prevent  range  smearing  the  following 
equation  must  be  satisfied 


2 


< 


(2.11) 


where  v  is  the  target  velocity  where  the  <  is  interpreted  as  less  than  by  a  factor  of  at 
least  10.  Solving  Equation  (2.11)  for  the  velocity  in  terms  of  the  time-bandwidth 
product  yields 
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V  <  _  (2.12) 

BT, 

r  d 

Therefore,  the  maximum  time-bandwidth  product  of  a  signal  is  limited  by  the 
maximum  target  velocity  of  interest. 

According  to  Equation  (2.12),  if  a  radar  is  to  process  targets  that  have  a 
maximum  velocity  of  v^,  then  the  time-bandwidth  product  should  be  chosen  to  be 

5  r  «  .1  —  (2.13) 

rd 

max 

Substituting  this  value  of  into  Equation  (2.6)  yields 


As  seen  in  Equation  (2.14),  the  larger  the  velocity  of  interest,  the  lower  the  radar- 
interceptor  advantage.  Notice  that  the  maximum  achievable  time-bandwidth  product 
for  a  300  m/s  target  is  100,000.  This  is  equivalent  to  a  50  dB  processing  gain.  Also 
notice  that  the  maximum  possible  processing  gain  decreases  at  10  dB/decade  with 
increasing  velocity.  This  results  in  a  5  dB/decade  usable  range  improvement. 


Range-Doppler  Trade-Off 


It  is  of  interest  to  investigate  the  product  of  range  resolution  and  Doppler 


resolution.  The  range-Doppler  resolution  product  is  given  by 
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Av  M  = _ — _  (2.15) 

Since  the  speed  of  light  is  constant,  and  it  is  assumed  that  the  carrier  frequency  is 
fixed  by  the  radar  hardware,  the  range-Doppler  resolution  product  is  fixed  by  the  time- 
bandwidth  product.  Therefore,  for  a  fixed  time-bandwidth  product  range 
resolution  must  be  traded  for  Doppler  resolution  or  vice  versa.  So,  for  a  fixed  time- 
bandwidth  product,  a  high  range  resolution  radar  signal  will  have  to  come  at  the 
expense  of  reduced  Doppler  resolution. 

2.7.  The  Cross-Ambiguity  Function  and  Clutter  Analysis 
When  a  radar  transmits  a  signal,  the  reflections  do  not  come  from  the  target 
alone.  There  are  other  competing  reflections  that  tend  to  obscure  the  wanted  target 
reflection.  These  unwanted  reflections  are  typically  referred  to  as  clutter  returns. 

A  general  analysis  tool  that  is  used  to  describe  the  response  of  a  radar  to  any 
given  clutter  distribution  is  the  cross-ambiguity  function  [3]  of  the  radar.  The  cross¬ 
ambiguity  function  describes  the  output  of  a  radar  to  a  point  target  located  at  any 
arbitrary  range  delay  x  and  Doppler  frequency  fp  and  is  given  by 

00 

X^('c  ./d)  =  ^  [  x{t)  y(r-T)  dt  .  (2-16) 

Properties  of  the  cross-ambiguity  function  will  be  discussed  in  detail  in  Chapter  HI. 

The  ideal  radar  signal  would  have  an  impulsive  cross-ambiguity  function.  In 
other  words,  targets  at  the  appropriate  range  and  velocity  would  give  a  maximum 
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response,  while  aU  other  targets  would  yield  no  response.  The  ideal  signal,  however, 
is  impossible  to  generate.  Stewart  and  Westerfield  [4]  investigated  the  cross¬ 
ambiguity  functions  four  general  signal  classes:  long  pulses,  short  pulses,  linear 
frequency  modulated  pulses,  and  pseudo-random  coded  pulses.  Density  plots  of  these 
four  signal  cross-ambiguity  functions  are  sketched  in  Figure  2.6.  The  shaded  portions 
of  the  plots  represent  the  range-Doppler  regions  to  which  the  radar  will  respond.  The 
long  pulse,  short  pulse,  and  frequency  modulated  cross-ambiguity  densities  are 
confined  to  a  single  ellipsoidal  mainlobe  that  is  l/B^  in  width  along  the  range  axis  and 
UTa  in  width  along  the  Doppler  axis.  The  ellipsoid,  however,  spans  over  a  portion  of 
the  surface  wide  in  Doppler  and  wide  in  range.  It  should  be  noted  that  for  the 
uncoded  signals  (the  long  pulse  and  short  pulse)  that  =  1  (i.e.,  no  pulse 
compression).  The  pseudo-random  cross-ambiguity  function  density  has  a  central 
spike  with  a  width  of  l/B^  along  the  range  axis  and  l/T^  along  the  Doppler  axis. 
Additionally,  it  has  a  uniform  pedestal  region  that  is  l/B^^  high  and  spans  over  a 
width  of  Br  in  Doppler  and  in  range.  In  other  words,  the  long  pulse,  short  pulse, 
and  frequency  modulated  radars  will  respond  strongly  to  targets  or  clutter  at  range  and 
Doppler  values  within  the  mainlobe  only.  For  the  pseudo-random  signal,  the  radar 
will  respond  strongly  to  targets  at  the  appropriate  mainlobe  range  and  Doppler,  and 
wUl  respond  in  a  weak  random  manner  for  targets  or  clutter  in  the  pedestal  region. 

Typically,  it  is  assumed  that  the  wanted  target  is  a  point  at  the  origin  of  the 
range-Doppler  plane,  while  clutter  is  distributed  in  range  and  Doppler.  Most  often,  it 
is  assumed  that  clutter  is  localized  in  Doppler  and  distributed  along  the  range  axis  as 


Figure  2.6.  Generic  Cross-Ambiguity  Functions  for  Four  Signal  Types. 


shown  in  Figure  2.7.  A  good  measure  of  clutter  rejection  for  this  clutter  distribution 
(for  waveform  comparison  purposes  only)  is  therefore,  the  integrated  sidelobe  level 
(ISL) 


23 


ISL  (4)  -  f  I  .  Q  p  *  ■ 

-00 

The  integrated  sidelobe  level  is  the  clutter  rejection  figure  of  merit  that  will  be  used 
throughout  this  study.  It  is  not  a  generalized  figure  of  merit,  and  does  not  apply  to 
every  specific  clutter  distribution.  It  does,  however,  serve  as  a  useful  figure  of  merit 
for  the  purposes  of  comparing  the  performance  of  different  waveforms. 

2.8.  Radar  Scattering  Phenomenology 

The  scattering  of  radar  waves  from  targets  is  characterized  by  three  different 
regions  that  are  dependent  upon  the  ratio  of  the  target’s  size  L  to  the  radar  wavelength 
X.  The  Rayleigh  region  is  associated  with  target  lengths  that  are  much  less  than  a 
wavelength,  or  L  <X.  The  resonance  region  is  associated  with  targets  whose  lengths 
are  on  the  order  of  a  wavelength  L~X.  The  optical  region  is  associated  with  targets 
that  are  much  larger  than  a  wavelength  L>X. 

2.8.A.  Rayleigh  Region  Scattering  Characteristics 
In  the  Raleigh  region,  scattering  phenomena  are  described  by  the  quasi-static 
dipole  moment  of  the  target  [5].  In  general,  the  scattered  field  changes  with  square  of 
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Figure  2.7. 
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Clutter  Distribution  Along  a  Constant-Doppler  Line. 
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the  frequency  of  the  excitation 


E  fo  (2-18) 

Since  the  wavelength  is  much  larger  than  the  target  length,  Rayleigh  scattered  returns 
have  the  advantage  of  being  aspect  independent  but  they  convey  little  target-specific 
information.  Additionally,  low  frequency  signals  make  it  difficult  to  perform  Doppler 
processing  and  high-gain  narrow-beamed  antennas  are  almost  impossible  to  build  since 
the  half-power  beamwidth  is  given  roughly  as  the  ratio  of  wavelength  to  aperture 
length  (Equation  (2.2)).  As  an  example  of  a  Rayleigh  scattering  target,  a  target  that  is 
15  meters  in  length  would  exhibit  Rayleigh  scattering  for  firequencies  lower  than  2 
MHz.  Therefore,  the  Rayleigh  region  is  not  of  interest  for  target  identification 
purposes. 


2.8.B.  Resonance  Region  Scattering  Characteristics 
The  resonance  region  is  defined  by  the  ratio  of  target  length  to  wavelength  on 
the  order  of  1  so  that 


1  <  £  <  10  (2.19) 

X 

Resonant  region  scattering  is  typically  described  by  highly  damped  resonant  modes. 
These  modes  result  from  the  waves  circumnavigating  the  various  portions  of  the  target. 
Therefore,  the  frequencies  of  the  modes  are  related  to  the  physical  size  of  the  different 
scattering  structures.  Theoretically,  since  the  resonant  modes  are  related  to  the 
physical  scattering  structure  dimensions,  resonant  mode  scattering  is  excellent  for 
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target  identification  purposes  [6,  7,  8,  9,  10].  The  mode  locations  are  relatively  aspect 
dependent  since  the  scattering  structures  are  still  a  small  number  of  wavelengths  in 
physical  dimension.  A  radar  to  excite  resonant  modes,  however,  would  have  similar 
problems  to  a  Rayleigh  region  radar.  Since  the  wavelength  is  on  the  order  of  the 
target’s  physical  dimension,  the  radar  antenna  would  have  to  be  much  greater  than  the 
target  dimensions  for  the  antenna  to  achieve  a  narrow  beamwidth  and  high  gain. 
Additionally,  Doppler  resolution  would  not  be  sufficient  for  clutter  rejection.  As  an 
example  of  a  resonant  mode  radar,  the  radar  would  operate  in  a  band  from  20  MHz  to 
200  MHz  for  a  15  meter  target.  From  Equation  (2.2),  the  radar  antenna  would  have  to 
be  on  the  order  of  200  meters  per  dimension  to  achieve  a  30  dBi  gain  at  the  low  end 
of  the  frequency  spectrum. 

2.8.C.  Optical  Region  Scattering  Characteristics 
The  optical  region  is  defined  by 

il  >  1  (2.20) 

A 

Optical  region  scattering  is  defined  by  impulses  in  the  time  domain  and  complex 
exponentials  in  the  frequency  domain.  The  location  of  the  impulses  in  the  time 
domain  are  associated  with  physical  scattering  centers  of  the  target.  These  scattering 
centers  vary  rapidly  with  target  aspect  angle,  but  they  do  provide  good  features  for 
target  identification  purposes  [11,  12].  Since  the  wavelength  is  small,  radar  antenna 
sizes  can  be  reduced  and  Doppler  processing  can  be  used  for  clutter  removal.  The 
disadvantage  is  that  the  scattering  centers  are  highly  aspect  dependent  (more  so  at 
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higher  and  higher  frequency),  and  therefore,  a  large  amount  of  data  is  required  to 
describe  the  target.  The  aspect  dependence  is  given  by  the  angular  sampling  criterion 
[11] 


AS  <  A  .  (2.21) 

4L 

Equation  (2.21)  is  required  for  coherent  data  processing.  For  noncoherent  processing, 
the  criterion  can  be  relaxed  to  [11] 

A0  <  -^  (2.22) 

L 

where  Ai?  is  the  range  resolution  of  the  radar.  Coherent  scattering  data  provides  a 
better  feature  vector  for  target  identification  than  noncoherent  data  [12]  but  at  the 
expense  of  a  much  larger  data  storage  requirement.  As  an  example,  consider  a  15 
meter  target  and  a  .15  meter  wavelength  (corresponding  to/p  =  2.0  GHz).  The 
angular  sampling  criterion  of  Equation  (2.21)  says  that  for  the  target  to  be  sampled 
properly  in  angle,  the  target  signature  must  be  sampled  every 

A0  < _ _  =  0.0025  radians 

(4)  (15)  (2.23) 

=  0.14° 

Therefore,  the  high  antenna  gain  and  Doppler  processing  capability  come  at  the 
expense  of  a  rapidly  varying  target.  For  a  15  meter  target  and  1  meter  resolution, 
noncoherent  processing  would  require 
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A0  <  _L  =  0.067  radians  ^2  24) 

‘i  3.84»  . 

Therefore,  coherent  processing  means  that  3.84/0.14  =  27.4  more  data  is  required  for 
database  storage  purposes  per  angular  dimension  for  a  target  identification  algorithm. 
To  describe  the  entire  target,  (27.4)^  =  751  times  more  data  is  required  for  coherent 
data.  In  general,  the  ratio  of  coherent  discritization  to  noncoherent  discritization  is 
given  by 

^  ^  (2.25) 

AS.  X 

Therefore,  the  larger  the  number  of  wavelengths  per  range  resolution  cell,  the  more 
data  will  be  required  for  target  storage. 

It  is  important  to  note  that  target  identification  requires  a  signal  bandwidth 
sufficient  to  resolve  individual  scattering  centers  on  a  given  target  (i.e.,  UWB 
waveforms).  This  typically  has  not  been  assumed  in  studies  aimed  at  improving  signal 
to  clutter  ratios.  For  this  reason,  signal  to  clutter  ratios  have  been  improved  only  by 
suppressing  clutter  returns  (lowering  the  ISL).  It  seems  that  a  logical  fusing  of  signal 
to  clutter  enhancement  and  target  identification  should  include  enhancing  the  target 
return  as  well  as  suppressing  the  clutter  return.  This  enhancement  could  be 
accomplished  by  designing  signals  with  cross-ambiguity  functions  with  mainlobes  that 
cause  the  scattering  centers  to  sum  coherently  rather  than  simply  measuring  the 
location  and  energy  from  each.  This  enhancement  would  have  to  come  without 
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increasing  the  ISL  in  a  like  manner.  This  concept  will  be  addressed  in  detail  in  this 
study. 

2.9.  Conclusions 

The  range  at  which  an  LPI  radar  obtains  an  advantage  over  an  intercept 
receiver  is  extended  by  increasing  the  time-bandwidth  product  of  the  radar  waveform 
and  by  enhancing  the  reflectivity  of  the  target.  The  time-bandwidth  product  is  limited 
by  maximum  velocity  of  interest  to  the  radar,  therefore  the  radar-interceptor  advantage 
gained  from  increasing  the  time-bandwidth  product  is  limited  by  target  velocity.  The 
next  logical  step  in  improving  the  radar-interceptor  advantage  is  therefore,  target 
enhancement  through  the  use  of  target-specific  information. 

The  clutter  rejection  performance  of  a  waveform  is  quantified,  for  general 
comparative  purposes,  by  the  integrated  sidelobe  level  (ISL)  of  the  cross-ambiguity 
function.  This  is  true  specifically  for  the  case  of  range-distributed  Doppler¬ 
concentrated  clutter.  Additionally,  the  cross-ambiguity  of  a  pseudo-random  LPI 
waveform  is  comprised  of  an  impulse-like  target-specific  and  a  uniformly  distributed 
pedestal.  The  pedestal  region  is  the  primary  source  of  clutter  contamination. 

Finally,  to  implement  a  realistic  radar  target  identification  technique,  the  only 
practical  target  scattering  regime  is  the  optical  region  because  of  the  requirements  on 
Doppler  resolution  and  antenna  gain.  This  advantage  comes  at  the  expense  of  an 
increasingly  complex  target  scattering  description.  Additionally,  the  target  description 
for  a  coherent  processor  is  much  more  complex  that  of  a  noncoherent  processor 
assuming  that  the  range  resolution  of  the  radar  is  many  wavelengths. 
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CHAPTER  m 


ANALYSIS  OF  A  GENERAL  FSK/PSK  RADAR  SYSTEM 

In  this  chapter,  the  radar  cross-ambiguity  function  [1,  p.  153-157]  wiU  be 
discussed.  The  cross-ambiguity  function  describes  the  response  of  a  radar  system  to 
an  impulse-like  target  (or  point  target)  located  at  an  arbitrary  range  and  Doppler  shift. 
In  this  sense,  the  cross -ambiguity  function  can  be  thought  of  as  the  impulse  response 
of  the  radar.  The  cross-ambiguity  function  is  the  general  case  of  the  radar  ambiguity 
function  [1,  pp.  118].  The  ambiguity  function  is  also  referred  to  as  the  matched-filter 
response,  and  the  uncertainty  function  [2].  The  cross-ambiguity  function  is  also 
related  to  the  cyclic  cross-correlation  function  as  discussed  in  [3,  pp.  369-373]. 

In  Section  3.1,  general  properties  of  the  cross-ambiguity  function  that  are 
particularly  applicable  to  radar  system  analysis  is  be  derived  First,  the  cross¬ 
ambiguity  function  of  two  signals  is  defined  along  with  a  discussion  of  the  underlying 
assumptions  made  during  its  derivation.  Second,  the  volumetric  distribution  of  the 
cross-ambiguity  function  is  analyzed  for  a  general  system.  Third,  the  property  of 
linearity  is  proven.  Fourth,  the  effect  of  convolving  one  of  the  two  signals  with  an 
impulse  response  on  the  cross-ambiguity  function  is  investigated.  Fifth,  the  effect  of 
signal  repetition,  or  periodicity,  on  the  cross-ambiguity  function  is  studied.  Finally, 
the  effect  of  additive  stationary  noise  on  the  cross-ambiguity  function  is  determined. 
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In  Section  3.2,  a  general  FSK/PSK  radar  system  diagram  is  presented  that 
includes  the  effects  of  the  target  transfer  function,  delay,  and  motion.  This  system 
diagram  is  analyzed  using  the  properties  derived  in  Section  3.1.  In  conclusion,  is  be 
seen  that  the  interaction  between  the  target  and  the  radar  is  defined  by  the  correlation 
of  the  radar  cross-ambiguity  function  with  the  down-converted  target  transfer  fimction. 


3.1.  The  Cross- Ambiguity  Function 


In  this  section,  the  properties  of  the  cross-ambiguity  function  are  discussed. 

The  cross-ambiguity  function  of  a  radar  is  a  rigorous  mathematical  description  of  a 
radar’s  response  to  an  ideal  point  target  moving  at  a  constant  range  rate.  The  cross¬ 
ambiguity  function  is  therefore  a  two  dimensional  fimction  of  range  delay  x  and 
Doppler  frequency /p.  The  cross-ambiguity  function  Xxyi'^fo)  the  signal  x(t)  with  the 
signal  y(t)  is  defined  as 


00 

.  /d  )  =  ^  [  yit-x)  dt  (3-1) 

where  is  the  duration  of  the  signal  x(t),  x  is  the  time  delay  between  waveforms,  and 
/o  is  the  Doppler  frequency  shift  induced  by  the  moving  target.  Equation  (3.1) 
describes  the  output  of  the  radar  receiver  for  various  values  of  t  and  Doppler 
frequency  Notice  that  Equation  (3.1)  is  identical  to  the  time  average  cross- 
correlation  function  [4]  except  for  the  Fourier  transform  kernel 


33 


3. LA.  Physical  Limitations  on  Resolution 
It  should  be  noted  that  the  cross-ambiguity  function  is  an  accurate 
representation  of  the  radar  response  under  certain  restrictions  on  the  signal  duration 
and  bandwidth  in  relation  to  the  target’s  velocity  v  and  acceleration  a.  These 
restrictions  are  derived  in  [1,  p.  61]  and  are  given  by 

(3.2) 

"  V 

and 

BTa  ^  0.2  £  (3.3) 

'  a 

where  c  is  the  speed  of  light  in  vacuum.  Therefore,  the  analysis  that  follows  is  only 
valid  within  the  restrictions  of  Equations  (3.2)  and  (3.3). 

An  example  calculation  will  now  be  performed  in  order  to  become  acquainted 
with  typical  velocity  and  acceleration  requirements  for  a  given  time-bandwidth 
product.  Consider  a  radar  with  a  bandwidth  of  B^  =  300  MHz  and  a  processing  time 
Tj  of  100  psec.  According  to  Equation  (3.2),  the  maximum  allowable  velocity  of  a 
target  is 

V  ^  0.1 - ^  ^  -  =  1000  m/sec  (3.4) 

(  300  X  10*"  )  (  100  X  10-®  ) 

and  the  maximum  allowable  acceleration  is 
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a  ^0  2 _ ^  ^  _  =  20  X  10"  m/sec^  .  (3.5) 

(  300  ;c  10"  )  (  100  .r  10-"  f 

The  acceleration  is  unreasonably  high  for  any  realistic  target  and  is  not  a  limiting 
factor,  but  the  velocity  limit  is  of  interest.  The  radar  signal  described  above  would  be 
limited  to  processing  reflections  from  targets  moving  at  velocities  of  less  than  1000 
m/sec  or  1943  kts.  It  should  be  noted  that  if  the  processing  time  is  reduced,  the 
maximum  target  velocity  will  increase. 

The  velocity  resolution  of  a  pulsed  radar  is  given  by 

Av  =  (3.6) 

J  o  a 

where /o  is  the  carrier  frequency  of  the  waveform.  For/^  =  2.0  GHz  and  =  100 
psec,  the  velocity  resolution  is 

Av  = _ ^  ^  _  =  750  m/sec  (3.7) 

2  (2.0  jc  10’)  (100  X  10-") 

which  is  on  the  order  of  the  maximum  velocity.  Therefore,  the  waveform  above 
cannot  resolve  the  target  velocity.  Notice,  however,  that  if  is  increased,  then  the 
velocity  resolution  is  increased. 

The  signal  described  above  is  well-suited  for  measuring  range  since  the  range 


resolution  of  a  signal  is  given  by 
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AR  =—  (3.8) 

25, 

again  where  B,  is  the  bandwidth  of  the  signal  (for  the  signal  bandwidth  of  above,  the 
range  resolution  is  1/2  meter).  The  product  of  the  range  resolution  and  Doppler 
resolution  is  given  by 

ARAv  = - — -  .  (3.9) 

4  /„  BT. 

•^0  r  a 

Therefore,  for  a  fixed  time- bandwidth  product  5,r^  and  carrier  frequency  /„,  range 
resolution  must  come  at  the  expense  of  Doppler  resolution  and  vice  versa.  Therefore, 
a  radar  with  high  range  resolution  cannot  also  have  high  Doppler  resolution. 

Several  important  fundamental  properties  of  the  cross-ambiguity  function  are 
now  discussed. 


3.I.B.  Volume  Distribution  of  the  Cross- Ambiguity  Function 
It  is  of  interest  to  determine  the  volume  beneath  the  surface  of  the  ambiguity 
function.  This  volume  is  defined  by 


V  = 


J  J  i  fo  ) 


df^ 


(3.10) 


Following  the  method  shown  in  [1,  pp.  120-122],  the  volume  can  be  shown  to  equal 
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y  ■  xjo,o)  Xy,m  (311) 

where 

Xj0,0)  =±  (3-12) 

^  d  ^ 

is  the  normalized  energy  in  x(0  and 

x„(o,o)  =  ±  r  I  Ko  r  lA 

^  d  ^ 

is  the  normalized  energy  in  y(t).  Therefore,  for  fixed  signal  power,  the  volume 
beneath  the  cross-ambiguity  function  is  fixed.  Suppressing  the  surface  of  the  cross¬ 
ambiguity  function  in  one  particular  region  results  in  an  increase  elsewhere  on  the 
surface. 

Two  more  properties  of  the  cross-ambiguity  function  are  related  to  the 
distribution  of  the  volume  underneath  the  function.  The  first  relates  the  delay-axis 
function  to  the  volume  distribution  in  Doppler: 


)l^  dx  =  r  xLC-t.O)  x^(x,0)  dx  .  (3-14) 

-00 

The  left  side  of  Equation  (3.14)  is  the  integrated  sidelobe  level  (ISL).  Therefore,  the 
ISL  is  the  Fourier  transform  of  the  product  of  the  autocorrelation  functions  of  the  two 
signals  used  in  the  signal  processor.  From  the  Wiener-Khinchin  theorem,  the 
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integrated  sidelobe  level  is  equal  to  the  convolution  of  the  two  signal  power  spectral 
densities  such  that 

ISL  ( /,  )  "  S^(  /,  )  *  y  /,  )  (3.15) 

Therefore,  the  ISL  of  a  radar  signal  processing  scheme  is  dependent  more  upon  the 
power  spectral  densities  of  the  signals  used  in  the  processor  than  the  actual  signal  type 
(e.g.,  linear  frequency  modulation  versus  pseudorandom  modulation). 

The  second  property  relates  the  Doppler-axis  function  to  the  volume 
distribution  in  delay; 


00  00 

4 ')  I "  =  f  xi(0,  fo )  yo,  /„ )  .  (3-16) 

Equation  (3.16)  shows  that  a  signal  with  a  narrow  range  mainlobe  will  have  energy 
spread  over  a  broad  region  in  Doppler  (broad  bandwidth),  and  likewise  a  signal  with  a 
broad  range  mainlobe  will  have  its  Doppler  energy  concentrated  in  a  narrow  Doppler 
mainlobe  (narrow  bandwidth).  It  should  be  noted  that  the  majority  of  the  volume 
beneath  the  cross-ambiguity  function  of  a  pseudorandom  signal  is  contained  beneath 
the  pedestal  region. 


3.I.C.  Linearity 

The  property  of  linearity  is  one  of  the  most  important  properties  in  systems 
analysis.  This  property  makes  it  possible  to  easily  analyze  systems  with  multiple 
inputs  such  as  target  return,  clutter,  and  noise.  The  cross-ambiguity  function  of  t) 
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and  y(t)  is  a  linear  operation  on  x(t)  and  y(t).  This  can  be  seen  by  setting  y(t)  =  y,(t) 
+  y2(t)  which  results  in 


00 

.  4  )  =  ^  f  ^(0  yit-T)  dt 
=  X:^,(^  .4  )  ^  X;.,/  "C  .4  ) 


(3.17) 


3.I.D.  Linear  System  Effects 

A  particularly  useful  property  of  the  cross- ambiguity  function  is  the  case  for 
which  the  signal  y(t)  has  been  convolved  with  a  filter  impulse  response  h(t) 


X;^.fc(''’/D  )  =  -^  f  xiO  (yit-x)  *h{t))  dt  .  (3-18) 

This  property  makes  it  possible  to  model  distributed  targets  with  a  general  impulse 
response  function  h(t).  This  impulse  response  h(t)  is  referred  to  as  the  range  profile  of 
the  target  since  it  describes  how  the  target  behaves  in  the  time  (range)  domain. 
Applying  the  convolution  operation  to  Equation  (3.18)  yields 


00 

f  ” 

r  x(o 

\yit-x-u)  h{u)  du 

do 

dt 


(3.19) 


After  switching  the  order  of  integration,  the  cross-ambiguity  function  becomes 
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00/00 

4  )  =  ^  I  J  40  3'(r-T-:^) 


</r 


So  I 


h(,u)  du 


(3.20) 


which  can  be  seen  to  equal 


4  )  =  ^  [  (  4  )  ) 

Therefore,  the  interaction  between  the  radar  and  a  target  is  described  completely  by 
the  cross-ambiguity  function  and  the  target  impulse  response. 


3.I.E.  Waveform  Repetition 

Waveform  repetition  is  often  used  in  radar  systems  to  create  a  clutter-clear  area 
around  the  mainlobe  of  the  cross-ambiguity  function  mainlobe  along  the  Doppler-axis 
[1,  pp.  141-143].  Signal  periodicity  can  be  modeled  as  a  convolution  of  a  short- 
duration  repeated  signal  x/t)  with  a  periodic  train  of  impulses.  For  an  ideally  periodic 
train  of  impulses,  there  would  be  infinitely  many  impulses.  But  for  a  study  of  realistic 
signals,  a  finite  duration  impulse  train  is  used.  A  periodic  signal  can  therefore  be 
written  as 


40  ■'”4^ 

m«0 


(3.22) 


where  Tp  is  the  period  of  the  signal  and  M  is  the  number  of  repetitions  of  the  short¬ 
term  waveform.  Equation  (3.22)  can  also  be  written  in  terms  of  a  periodic  impulse 
response  function  hp( t)  as 
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xit)=x^{t)*h^it)  (3.23) 

Writing  the  cross-ambiguity  function  in  terms  of  the  periodic  impulse  response 
fimction  hp(t)  yields 


OO  00  00 


^  d  -00  ^  do 

hpiu^)  hpi,)  du^  du^  . 


(3.24) 


After  a  change  of  variable  of  z=t-Uj,  Equation  (3.24)  becomes 


00  00  CO 


X^(^,  4  )  =  I  J  J  ^.(2)  }',(^-'c-(«2-Wj))  dz 


-00  -CO  -CO 


(3.25) 


^(“i)  ^(“2)  <^“1  ^“2  • 

where  the  identity  MTp  =  has  been  used.  By  making  the  substitution  of  the  cross¬ 
ambiguity  function  of  the  short  term  waveform  imbedded  within  the  integral  of 
Equation  (3.25)  yields 


00  00 

X;^('^.  4  )  =  ^  J  J  X.^,(-^-(M2-«i)  .  4  ) 

Ho  -CO 

/i  («j)  h  iu^  JUj  du^  . 


(3.26) 


Performing  the  convolutions  and  making  use  of  the  sifting  property  of  Dirac  delta 
functions  yields  the  following  summation 
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x„(t,  /d  )  =  4  S  2  ,  /„  )  e 

m«0  n=0 


(3.27) 


Using  the  substitution  of  I  =  m  -  n  casts  Equation  (3.27)  into  the  separable  summation 


M-l  M-l 

4  )  =  -4  E  X.^/^+^7’,  .  4  )  E  ^ 

M  /aO  m«; 


(3.28) 


-1 


M-l+l 


*  E  x,,,(t-'7-,  ,  /„ )  e  ^ 

^  (Af-1) 


/2n4m7' 


m*0 


The  summations  over  m  can  be  written  in  closed  form  by  using  known  properties  of 


geometric  series 


^  gnrfjnT  _ 
m»0 


(3.29) 


1 

and  Euler’s  identity.  Therefore,  Equation  (3.28)  can  be  written  as 


x>/) =4  ^ 


sin(7T/  {M  \l\)T^) 
sin{nf  rp 


(3.30) 


Note  that  the  /  index  locates  the  periodicities  of  the  mainlobes  of  the  cross-ambiguity 
function  of  the  short-term  functions  xjt)  and  yjt). 


3.I.F.  Noise  Contamination 

The  effect  of  additive  stationary  noise  is  now  quantified.  Since  the  cross¬ 
ambiguity  function  is  a  linear  transform,  the  noise  can  be  dealt  with  separately  from 


the  desired  receive  signal.  Therefore,  the  noise  component  of  the  received  signal  is 
given  by 
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^  f  nit-x)  dt  (3-31) 

da 

where  n(t)  is  a  stationary  noise  contaminant.  Assuming  that  the  noise  contaminant  is  a 
zero  mean  process,  the  expected  value  of  the  cross -ambiguity  function  between  the 
reference  signal  and  the  noise  contaminant  is  zero: 


4  )  }  =  ^  r  40  £  {  .  (3.3: 

^d  da 

A  more  important  measure  of  the  effect  of  noise  on  the  cross-ambiguity 
function  is  the  autocorrelation  fimction  of  the  noise  contamination.  This  is  given  by 


^  4  )  4)  }  =  A  [  [ 

T/  ll  (3.33) 

E  {  n(t-Xj)  n\t'-x^  } 

where  n{t}  is  a  stationary  noise  process.  By  making  the  substitution  of  t=z+t’  in 
Equation  (3.33),  the  following  equation  results 
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E  {  xjt,.  4  )  X„(v  /d  )  }  '  ^  J  J  j:  ■(( ')  41 ' «)  * ' 

7^  ico  -00  (3.34-} 

where  i?„„('zj  is  the  autocorrelation  function  of  the  noise  contaminant.  Noting  the 
autocorrelation  function  of  xftj  in  Equation  (3.34)  and  making  the  substitution  of 
A=T2-Tj,  the  autocorrelation  function  of  the  cross-ambiguity  function  equals 

00 

E  {  X..(X,  /„  )  4  )}  =  ^  f  EJZ)  Jtjz-A)  e'^'  dz  .  (3.35) 

Therefore,  the  autocorrelation  function  of  the  noise  takes  the  shape  of  the  convolution 
of  the  signal  autocorrelation  function  with  the  autocorrelation  function  of  the  noise. 
Equation  (3.35)  will  be  evaluated  for  two  special  cases  that  are  of  particular  interest: 
stationary  white  Gaussian  noise  and  stationary  narrowband  noise.  The  white  noise 
process  is  typical  of  thermal  noise  which  might  occur  in  the  radar  receiver  while  the 
narrowband  noise  is  typical  of  an  interference  type  signal  such  as  an  in-band  emitter. 

White  noise  is  defined  by  an  autocorrelation  function  =  Ngb(z)  where  Ng 
is  the  noise  single  sided  power  spectral  density.  Under  the  specific  case  of  white 
noise,  the  autocorrelation  function  of  the  noise  contaminant  is  given  by 

E  {  4  )  4  )  }  -  ^  «„(A)  e  .  (3.36) 

^  d 

Therefore,  the  noise  autocorrelation  assumes  the  shape  of  the  reference  signal 
autocorrelation  function.  Notice  that  Equation  (3.36)  is  not  a  function  of  t.  This 
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means  that  the  noise  variance  is  spread  uniformly  over  the  entire  delay-frequency 
plane.  Notice  also  that  the  noise  variance  decreases  as  increases  for  a  constant 
mean  square  signal  value  RJO).  Therefore,  for  a  constant  envelope  signal,  the  noise 
variance  is  a  function  of  the  processing  time  and  additive  noise  variance  only. 

The  narrowband  noise  contaminant  case  is  defined  by  RJiz)  =  cos(2nf„  z). 
Narrowband  noise  manifests  itself  in  the  following  manner 

e  { /» )  /„  ) )  =  ^  ( sj  4  - 

*  SJ  ) 

where 

-00 

is  the  power  spectral  density  of  the  signal  x(t).  Notice  that  narrowband  noise  does 
vary  as  a  function  of  frequency.  It  causes  a  contamination  that  has  the  same  spectral 
shape  as  the  signal  power  spectral  density  except  for  a  simple  modulation  by  the  noise 
contamination  frequency.  In  the  delay  dimension,  there  is  a  phase  shift  only  for 
complex  signals.  For  a  real  signal  along  the  zero  Doppler  axis  the  shape  of  the  noise 
carrier  wave  is  seen. 

As  an  example  of  narrowband  noise  contamination,  consider  the  case  of 
narrowband  contaminant  with  and  with  the  signal  power  spectral  density  S„(f) 

uniformly  distributed  in  a  bandwidth  from  -BJ2  to  BJ2  and  has  a  uniform  amplitude 
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of  HB.  This  is  shown  graphically  in  Figure  3.1.  The  interference  variance  at/o=0  is 
equal  to 


vfl'-  {  X  J'c.O)  } 


BT, 


(3.39) 


Therefore,  FSK/PSK  signaling  spreads  the  narrowband  signal  over  a  bandwidth  of  , 
and  the  contamination  level  is  improved  by  a  factor  of  B^^  . 


3.2.  General  FSK/PSK  Radar  System  Diagram 
A  generalized  radar  block  diagram  is  shown  in  Figure  3.2.  This  system 
implements  the  cross-ambiguity  function  of  the  transmit  signal  y(t)  with  a  reference 
signal  x(0.  For  an  FSK/PSK  radar,  the  transmit  signal  is  defmed  by 


iV-I 


y(0  =  E  n 


t-mT 


m«0 


cos  (271  (4,+/,)r) 


(3.40) 


where  a„  is  either  a  i  or  is  the  baseband  frequency  of  the  m"'  sub-pulse,  andf,  is 
the  carrier  frequency.  The  reference  signal  is  defmed  by 


N-l 


=E 


t-nT 


R>0 


,  -;■(  271 


(3.41) 


where  a„  is  either  a  /  or  -/,/„  is  the  baseband  frequency  of  the  n"'  sub-pulse,^  is  the 
carrier  frequency,  is  the  Doppler  correction  frequency,  and  is  the  phase  of  the  n"‘ 
sub-pulse. 


In  general,  the  transmit  signal  is  delayed  t  seconds  because  of  the  propagation 
through  the  atmosphere  at  the  speed  of  light  c.  The  signal  is  also  convolved  with  the 
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range  profile  of  the  target  h( t)  and  Doppler  shifted  by  ^  Hz  because  of  the  target 
motion.  With  these  considerations  taken  into  account,  the  reflected  signal  at  the  input 
of  the  receiver  is  given  by 


rit)  = 

N-\ 

E 

n 

'  t-mT-x  ] 

T 

^  m*0 

1  ^  J 

(3.42) 


X  cos  (  2k  (/,*/,  )<  -  2jt  (/.-/,)  T  )).  /,(0  . 

Referring  back  to  Figure  3.2,  it  is  seen  that  the  output  of  the  radar  system  is  given  by 


00 

fo)  f 

^  d  ^ 


(3.43) 


Substituting  Equations  (3.41)  and  (3.42)  into  Equation  (3.43)  along  with  the  principles 
derived  in  Section  3.1  results  in 


/d  )  =  f  (  fo  -/c  )  )  Ku)  du 

J 


(3.44) 


where 


N-l  N-l 


/d  -  /c  )  =  E  E  X„„('C.  fo  -fc) 

m»0  n*0 


(3.45) 


The  function  within  the  summation  represents  the  cross-ambiguity  function  of  each  of 
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the  frequency  sub-pulses  of  the  signal  and  is  equal  to 


,  .  f  .  1  r  TT  f  yj  ( t-nT 

fo  -fc)  =  11  - f—  11  -y- 

dc  [  ^  J  [ 


(3.46) 


Evaluating  the  integral  of  Equation  (3.46)  yields 


Xm„(  irn-n^E)  T  ,  ^  (1-lel) 


X  sinc[n(f^-f^^(f^-f„)  (1-lel)  r  } 

where  the  substitution  of  t  =  (m-n+e.)  T  has  been  made  and  -1  <  e  <  1.  For  lei  >  i, 
Equation  (3.47)  is  identically  equal  to  zero.  Making  use  of  this  representation  of  x, 
and  letting  I  =  (m-n).  Equation  (3.47)  can  be  expressed  as  the  following  single  index 


(3.47) 


summation 


Xxy(^+e,  fo  -fc)  = 


E  -/c  )  •"  E  f D  “ /c  )  (3.48) 

mW  mW+1  e![0,l] 

•  • 

AT+Z-l  W+/-2  T^O 

E  X„(„-o(e-/D  -fc)  ^  T,  X.(.-M)(-l-e>4  -/c  )  ^^0 

m’O  e:[-l,0] 

Notice  that  Equations  (3.46-3.48)  are  all  considered  to  be  baseband  terms  because  they 

are  all  dependent  on  the  baseband  frequency  terms  only.  Equation  (3.49)  shows  why 

this  is  true.  In  Equation  (3.49),  the  carrier  term  ^  is  factored  out.  Writing  Equation 

(3.49)  in  terms  of  the  correlation  integral  shows  that  this  separation  of  the  carrier  term 
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results  in  the  down-conversion  of  the  target  impulse  response  to  baseband: 

xA /„-/,)=  f  X;,(^*«. 4  - /J  ( KM  « )  ‘‘x 

-00 

3.3.  Conclusions 

The  volume  beneath  the  surface  of  the  cross-ambiguity  function  is  constant  for 
constant  power  signals.  Therefore,  reductions  in  ambiguity  function  volume  in  one 
region  result  in  increases  in  other  regions.  This  means  that  if  the  ISL  for  some 
specific  Doppler  frequency  is  reduced,  then  the  ISL  will  be  increased  for  other 
Doppler  frequencies.  Signal  repetition  causes  the  Doppler  response  of  the  repeated 
signal  to  be  multiplied  by  a  weighting  function  that  results  in  a  sharper  mainlobe  in 
Doppler  at  the  expense  of  range  ambiguities.  In  other  words,  the  ISL  is  reduced  in  the 
moving  target  region  of  the  cross-ambiguity  function,  but  the  stationary  target  ISL  is 
increased. 

The  response  of  a  radar  to  an  arbitrary  target  range  profile  is  given  by  the 
correlation  of  the  baseband  signal  cross-ambiguity  function  with  the  down-converted 
target  range  profile.  Therefore,  the  response  of  the  radar  to  a  given  target  is 
dependent  not  only  on  the  target  range  profile,  but  also  on  the  radar  signal  processor 
itself.  This  means  that  the  effective  radar  cross  section  of  a  target  can  be  enhanced  by 
designing  a  signal  processor  that  excites  an  maximum  response. 

Additive  white  noise  becomes  distributed  uniformly  over  the  entire  delay  and 
frequency  plane.  The  variance  of  the  additive  white  noise  decreases  as  the  processing 
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time  increases.  The  processing  time,  however,  can  only  be  increased  within  the 
confines  of  the  limitation  placed  on  the  time-bandwidth  product  by  maximum  target 
velocity.  Narrowband  noise  contamination  results  in  a  modulated  signal  spectrum 
contamination.  Therefore,  this  contamination  is  a  function  of  signal  spectral  shape.  If 
a  narrowband  contamination  is  introduced  to  the  receiver  signal  processor,  the  signal 
processor  tends  to  spread  the  energy  out  over  the  entire  bandwidth  of  the  signal 
(which  is  much  larger  than  the  contamination)  thus  reducing  the  effect  of  the  noise. 
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CHAPTER  IV 


PROBABILISTIC  ANALYSIS  OF  THE 
FSK/PSK  CROSS-AMBIGUITY  FUNCTION 

In  this  chapter,  the  FSK/PSK  cross-ambiguity  function  is  analyzed 
probabilistically  as  a  random  process.  The  function  is  analyzed  for  two  different 
regions;  the  mainlobe  region  and  the  sidelobe  (or  pedestal)  region.  The  mainlobe 
region  is  that  portion  of  Equation  (3.48)  for  which  1=0  and  IJ]  <  IIT^.  The  sidelobe 
region  is  that  portion  of  Equation  (3.48)  for  which  and/or  ///  >  IIT^.  Probabilistic 
analysis  of  these  functions  is  required  since  the  frequency  and  phase  sequences  within 
a  signal  are  inherently  random.  Therefore,  the  function  is  best  quantified  in  general  by 
its  mean  (or  expected  value)  and  variance.  The  expected  value  of  the  cross-ambiguity 
function  describes  the  shape  of  the  function  on  the  average  whereas  the  variance 
describes  how  the  realized  function  deviates  from  the  expected  value  on  the  average. 

It  is  assumed  that  the  frequency  and  phase  sequences  are  stationary  random  sequences. 
The  requirement  of  stationarity  is  necessary  for  the  analysis  that  follows  and  is  also  a 
valid  assumption  for  signals  that  are  random  from  signal-to-signal. 
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4.1.  Mainlobe  Analysis 

The  mainlobe  of  the  of  the  cross-ambiguity  function  of  an  FSK/PSK  radar 
signal  is  given  by 


N-l 


JV-l 


mM) 


m»l 


(4.1) 


mM) 


where  u(x)  is  the  unit  step  function  andf=fD-fc. 

Taking  the  expectation  of  the  sub-pulse  cross-ambiguity  functions  with  respect 
to  a„  and  a„  yields 


5mc{7t/(l  -lEl)r}  e  -^2,tr„er  (4.2) 

—  , 

0  rjv^n 

.. 

assuming  equal  likelihood  of  +1  or  -1  for  a„.  From  this  knowledge,  the  expectation  of 
Equation  (4.1)  with  respect  to  a„  and  a„  can  be  written  as 


N-l 


\a.{X^(^./)}  =E  X™„(^./) 


(4.3) 


m*0 


Now,  taking  the  expectation  with  respect  to  the  sub-pulse  frequencies  and  phases 
(])„,  the  expected  function  of  the  mainlobe  is  given  by 
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^  ^  sinc{nj{l -lel)r}  e 


ffi*0 


which  simplifies  to 


(4.4) 


(4.5) 

K^l 

smjnNjT)  ^jn(N-iw  p  g 
sin(n/r)  * 

where  K  is  the  number  of  oscillators  used  to  construct  the  signal,  is  the  probability 
that  the  k!'‘  oscillator  is  selected  within  the  burst,  and  (|)jt  is  the  phase  associated  with 
the  k!''  oscillator. 

The  mainlobe  variance  is  now  derived.  The  variance  is  important  because  it 
describes  the  mean  square  error  between  any  realization  of  an  FSK/PSK  cross¬ 
ambiguity  function  mainlobe  and  the  mean  mainlobe  function  given  in  Equation  (4.5). 
The  variance  of  the  mainlobe  of  the  cross-ambiguity  function  is  found  from  the 


following  equation 
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E  X„(„-i)(l-e./Me) 


m*l 


(4.6) 


+  E 


N-2 


E  X^„.n(-l-e./)  «(-e) 


mM) 


Substituting  in  the  individual  functions  into  Equation  (4.6)  and  collecting  terms  yields 
the  following 

{  X^(e.  /  )  }  = 

-L(l-lel)^  sinc^  {7i^l-lel)r} 

N 

+  £LR  lep  E  E  Pi  -f^  )  lei  T  } 

iM) 


2 

^r-1 

1  - 

Notice  that  the  variance  is  upper  bounded  by  1/N.  The  longer  the  sequence,  the  lower 
the  variance.  Therefore,  as  the  sequence  length  increases,  the  realized  cross-ambiguity 
function  converges  to  the  expected  function  Equation  (4.5). 
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4.2.  Sidelobe  Analysis 

The  sidelobe  region  of  the  FSK/PSK  radar  signal  cross-ambiguity  function  for 
positive  lag  values  is  given  by 


/ ) 


N-l 


E 


N-\ 


m»/+l 


(4.8) 


Using  Equation  (4.2),  the  expected  value  of  Equation  (4.8)  is  found  to  be  equal  to 
zero.  Therefore,  the  variance  of  Equation  (4.8)  is  given  by 


var  {  X^(/+e,  /  )  }  =  E  {  1  Xm(ni-/) 


mW+1 


Equation  (4.9)  can  also  be  expressed  as 


var{x^(/+e,/)  }  = 

(l-lel)^5^  Y,  PiP^sinc^ 

JW)  iM) 

^  E  E  PiPk ' { "(f-(/t-/;))iei^ }  • 

N  kM)  (^) 


(4.9) 


(4.10) 


Equation  (4.10)  was  derived  for  the  sidelobes  characterized  by  values  of  l>0.  It  is 
however  valid  for  both  positive  and  negative  values  of  /.  Again  notice  that  the 
variance  of  the  sidelobes  decreases  approximately  as  IIN. 
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4.3.  Integrated  Sidelobe  Level  Analysis 
A  more  important  result  obtained  from  Equation  (4.10)  is  the  integrate  sidelobe 
level  (ISL)  given  by 


ISL  1/T)  =  j  var  [  I  /  1  <  1/r  )  }  dx  .  (4-11) 


The  ISL  is  important  because  it  is  related  to  the  variance  (power)  of  the  ground  clutter 
return.  Substituting  Equation  (4.10)  into  Equation  (4.11)  and  evaluating  the  integral 
yields  the  following  equation  describing  the  ISL  in  terms  of  the  frequency  selection 
probabilities: 


/5L  (  I  /  I  <  1/r  )  =  2 


PjPk 


(4.12) 


Notice  that  the  ISL  is  not  a  function  of  the  sequence  length.  This  results  from  the  fact 
that  the  randomized  FSK/PSK  waveform  uniformly  spreads  its  sidelobe  energy  in 
range.  As  the  sequence  length  increases,  the  variance  of  the  sidelobes  decreases 
accordingly.  The  ISL  remains  constant  however,  since  the  domain  of  integration  is 
increased  in  a  proportionally. 

In  Chapter  HI,  the  effect  of  periodicity  (of  period  on  the  cross-ambiguity 
function  was  investigated.  From  the  results  of  Equation  (4.12)  and  (3.29),  it  can  be 
seen  that  the  ISL  of  an  FSK/PSK  waveform  that  is  repeated  M  times  is  given  by 


ISL^{\f  \<.IIT)  =  ISL  {f)  PER  {f  ,M)  (4.13) 


where 
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peru.m)’^  S 

M  (M-i) 


sin(K/(M  -  I/DTJ 
sin  (  n  /  r  ) 


(4.14) 


V  .  r  ,  J 

assuming  nonoverlapping  waveforms.  A  plot  of  Equation  (4.14)  for  M  =  2,  4,  8,  and 
16  is  shown  in  Figure  4.1.  As  is  shown  in  Figure  4.1,  the  zero  frequency  term  is 
increased  by  approximately  10  dB/decade  with  increasing  M  while  the  nonzero 
frequency  terms  are  decreased  approximately  as  10  dB/decade  with  increasing  M. 


4.4.  Conclusions 

The  mean  value  of  the  cross-ambiguity  function  is  non-zero  only  in  the 
mainlobe.  Additionally,  the  variance  of  the  cross-ambiguity  function  is  upper  bounded 
by  UN  at  the  mainlobe,  and  the  variance  falls  off  linearly  with  the  distance  from  the 
mainlobe  to  zero  at  its  outer  delay  boundary.  Even  though  the  variance  of  the 
sidelobes  can  be  decreased  by  increasing  the  sequence  length,  the  ISL  of  the  cross¬ 
ambiguity  function  is  unaffected  since  the  sidelobe  energy  is  simply  spread  over  a 
larger  area.  Additionally,  repetition  of  waveforms  has  the  effect  of  increasing  the  zero 
Doppler  region  ISL  whUe  the  moving  target  portion  is  decreased.  Therefore,  a  clear 
area  can  be  created  by  repeating  waveforms. 


Figure  4.1.  Periodicity  Function  in  dB  Versus  Frequency  for  Various  Values  of  M. 


CHAPTER  V 


MATCHED  FSK/PSK  SIGNAL  TIME-BANDWIDTH  PRODUCT 

In  Chapter  n,  it  was  seen  certain  factors  that  contribute  to  the  LPI  performance 
of  a  radar  signal  processor;  the  time-bandwidth  product  of  the  transmit  signal  and  the 
effective  radar  cross  section  of  the  target.  The  higher  the  time-bandwidth  product,  the 
better  the  LPI  performance  of  the  radar.  In  this  chapter,  the  effective  time-bandwidth 
product  of  an  FSK/PSK  radar  signal  is  derived.  The  time-bandwidth  product  is 
determined  by  first  calculating  the  effective  bandwidth  of  the  FSK/PSK  transmit  signal 
power  spectral  density.  The  effective  bandwidth  of  the  signal  power  spectral  density 
is  developed  based  on  the  signal’s  noise  equivalent  bandwidth  [1], 

5.1.  FSK/PSK  Radar  Signal  Power  Spectral  Density 
The  power  spectral  density  of  a  signal  describes  how  the  signal  power  is 
distributed  in  the  frequency  domain  on  the  average.  The  power  spectral  density  of  a 
signal  y(t)  is  defined  by  [1] 

S  («»  =  (5-1) 

^  d 

where  E{}  is  the  expected  value  operator,  is  the  signal  duration,  and  Tfo))  is  the 
Fourier  transform  of  y(t}  defined  by 
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7(co)  =  f  y{t)  e  dt  .  (5-2) 

-00 

Mathematically,  the  power  spectral  density  of  a  signal  is  equal  to  the  limit  of  Equation 
(5.1)  as  Tj  grows  to  infinity.  For  the  purposes  of  this  study,  however,  only  fmite- 
duration  signals  are  of  interest.  Therefore,  Equation  (5.1)  will  be  termed  the  finite- 
duration  power  spectral  density  of  the  transmit  signal. 

The  transmit  signal  of  the  radar  is  defined  by  Equation  (3.40).  Evaluating 
Equation  (5.3)  with  Equation  (3.40)  as  the  argument  of  the  integral  yields 
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Taking  the  expected  value  with  respect  to  the  random  FSK  code  and  the  random  PSK 
code  of  the  magnitude  squared  of  Equation  (5.3)  yields 


E  1  T(©) 


4  i-O 


sine" 


© 


-  (  )  J  ^ 


(5.4) 


NT 


K-l 


where  represent  the  probability  of  selecting  oscillator  k.  Now,  dividing  through 
by  the  duration  NT  of  the  signal  yields  the  finite-duration  power  spectral  density 
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Notice  that  the  power  spectral  density  of  the  signal  takes  the  shape  of  the  probability 
density  function  of  the  frequency  selection  process.  Mathematically,  the  term  fimte- 
duration  power  spectral  density  is  an  oxymoron.  By  definition,  power  signals  are 
infinite  in  duration  (persistent).  In  actuality,  y( t)  is  an  energy  signal  because  only  a 
finite  duration  portion  of  the  signal  is  processed  in  any  processing  interval.  For  the 
purposes  of  this  study,  it  will  be  assumed  that  the  actual  signal  exists  for  infinity,  but 
only  a  finite-duration  portion  of  the  signal  is  of  interest.  Therefore,  Equation  (5.5)  is 
valid  for  the  analysis  that  follows. 


5.2.  Equivalent  Bandwidth 

Equation  (5.5)  can  take  on  any  arbitrary  shape  that  has  the  properties  of  a  valid 
probability  density  function.  Therefore,  the  bandwidth  of  the  power  spectral  density  is 
a  difficult  parameter  to  define.  The  bandwidth  could  be  defined  to  be  (co/^.j-cog)  since 
that  is  the  maximum  possible  bandwidth  in  which  the  signal  could  possibly  exist.  A 
definition  of  filter  bandwidth  that  is  often  used  in  signal  processing  analysis  is  the 
noise  equivalent  bandwidth  of  a  signal.  A  pictorial  representation  of  noise  equivalent 
bandwidth  is  shown  in  Figure  5.1  and  mathematically  by  [1] 
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da 


(5.6) 
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Figure  5.1.  Equivalent  Bandwidth  of  a  Signal. 

Equation  (5.6)  assumes  that  the  oscillator  frequencies  are  separated  by  at  least  1/r  Hz. 
In  other  words,  the  different  frequency  sub-pulse  power  spectral  densities  do  not 
overlap.  From  Equation  (5.5),  the  maximum  value  of  the  power  spectral  density 
function  is  equal  to  one  foiu-th  of  T  times  the  maximum  probability  of  frequency 
occurrence 

max{  5^(co)  }  =  L  max{  }  .  (5.7) 

Additionally,  the  integral  of  Equation  (5.6)  is  in  a  form  that  is  tabulated  in  many 
integral  tables.  The  value  of  the  integral  can  be  foimd  to  equal 
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(5.8) 


Therefore,  according  to  Equation  (5.6),  the  equivalent  bandwidth  of  an  FSK/PSK  radar 
signal  is  given  by 


B  = 


(5.9) 


T  max|  } 

Equation  (5.9)  can  be  evaluated  for  two  extreme  probability  density  functions  in  order 
to  become  acquainted  with  the  equivalent  bandwidth.  The  first  density  function  is  a 
uniform  distribution  in  which  all  frequencies  occur  with  probability  IIK.  Under  this 
condition,  the  equivalent  bandwidth  is  given  by 


B  =  A  .  (5.10) 

.  r 

The  second  density  function  that  is  of  interest  is  one  in  which  only  a  single  frequency 
is  used.  Under  this  density  fimction,  the  equivalent  bandwidth  becomes 

S  .1  .  (5.11) 

'  T 

Both  of  these  equivalent  bandwidth  calculations  have  results  that  are  intuitively 
obvious. 

5.3.  Equivalent  Time-Bandwidth  Product 
Since  the  duration  of  the  FSK/PSK  radar  signal  is  equal  to  the  number  of 
pulses  time  the  duration  of  each  pulse,  the  equivalent  time-bandwidth  product  is  equal 
to 
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= 

e  a 


N 


(5.12) 


max{  p,  } 

It  should  be  noted  that  the  equivalent  time-bandwidth  product  should  be  used  to 
evaluate  the  LPI  performance  of  the  radar  signal.  It  should  not  be  used  to  determine 
the  maximum  target  velocity  that  can  be  processed  by  the  radar.  For  the  velocity 
limitation,  it  should  always  be  assumed  that  all  frequencies  will  be  utilized  equally 
often.  In  other  words,  for  determining  the  maximum  velocity  that  the  radar  can 
process,  the  following  equation  should  be  used 


BT=NK  .  (5.13) 

r  d 

It  should  also  be  noted  that  the  maximum  time-bandwidth  product  possible  occurs  for 
a  signal  that  uses  all  sub-pulse  frequencies  with  equal  likelihood.  Therefore,  as  a 
signal  takes  on  any  shape  other  than  a  uniform  shape,  the  equivalent  time-bandwidth 
product,  and  the  associated  LPI  performance  is  reduced.  A  measure  of  LPI 
performance  degradation  due  to  a  nonuniform  transmit  signal  power  spectral  density 
can  be  defined  as  the  ratio  of  Equation  (5.13)  to  (5.12) 


L 


LPI 


1 

K  max|  } 


(5.14) 


Notice  that  a  value  of  Lj^p,  less  than  unity  represents  an  LPI  performance  degradation. 
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5.4.  Conclusions 

The  power  spectral  density  of  an  FSK/PSK  radar  signal  is  defined  by  the 
probability  density  function  of  the  frequency  selection  random  process.  Since  the 
power  spectral  density  can  take  on  any  shape  that  is  a  valid  probability  density 
function,  the  bandwidth  of  the  signal  is  somewhat  difficult  to  parameterize.  For  the 
purposes  of  LPI  performance,  the  equivalent  bandwidth  of  the  radar  signal  is  used. 
This  equivalent  bandwidth  of  a  signal  is  equal  to  the  bandwidth  of  a  signal  with  equal 
power  to  the  original  signal  that  has  an  ideal  bandlimited  shape  and  a  magnitude  equal 
to  the  peak  value  of  the  actual  signal  power  spectral  density.  The  equivalent 
bandwidth  is  maximum  for  a  signal  that  utilizes  each  available  oscillator  with  equal 
likelihood  (uniform  power  spectral  density).  The  equivalent  bandwidth  is  bounded 
between  l/T  and  KIT  where  K  is  the  number  of  oscillators  available  and  T  is  the 
pulsewidth. 
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CHAPTER  VI 


MATCHED  FSK/PSK  SIGNAL  DESIGN 

In  this  chapter,  the  probabilistic  analysis  of  the  FSK/PSK  cross-ambiguity 
function  will  be  used  to  develop  design  techniques  for  Matched  FSK/PSK  radar 
signals.  Matched  FSK/PSK  radar  signals  are  FSK/PSK  signals  in  which  the  mainlobe 
of  the  cross-ambiguity  function  has  been  specifically  shaped  to  become  a  matched 
filter  for  a  given  target  range  profile.  The  design  technique  will  be  derived  from 
Equation  (4.5).  First  the  special  case  of  the  FSK/PSK  cross-ambiguity  function 
mainlobe  at  zero  frequency  will  be  analyzed.  This  wiU  be  shown  to  be  of  the  form  of 
a  discrete  Fourier  series.  The  discrete  Fourier  series  representation  wiU  allow  the 
mainlobe  of  the  FSK/PSK  cross-ambiguity  function  to  be  designed  to  approximate  the 
shape  of  any  arbitrary  time-limited  function.  A  step-by-step  example  involving 
computer  simulations  will  be  used  to  illustrate  the  concept. 

6.1.  FSK/PSK  Cross- Ambiguity  Function  as  a  Fourier  Series 
In  this  section,  the  mainlobe  of  the  FSK/PSK  cross-ambiguity  function  for  the 
case  of  zero  frequency  deviation  will  be  analyzed  to  show  that  it  is  of  the  form  of  a 
discrete  Fourier  series.  From  Equation  (4.5),  the  expected  value  of  the  mainlobe  of  an 
FSK/PSK  radar  signal  cross-ambiguity  function  is  given  by 
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E. /  )  }  -  ^  if/.c{7iAl-tel)7-)  e 

■  ■  (6.1: 

sin(K/V/7)  gjvJiN-i)fr  n  ^ 
sin(7t/7)  i=o  * 

again  where  K  is  the  number  of  oscillators  used  to  construct  the  signal,  p,.  is  the 
probability  that  the  /r""  oscillator  of  frequency  is  chosen  within  the  signal  burst,  and 
(j)jt  is  the  phase  of  the  /:'*  oscillator.  If  the  values  of  are  chosen  to  be  integer 
multiples  of  the  pulse  rate 


(6.2) 


then  the  summation  of  Equation  (6.1)  is  of  the  form  of  a  discrete  Fourier  series  (DFS). 
It  is  of  interest  to  express  Equation  (6.1)  in  a  more  convenient  form  of  the  product  of 
a  closed  form  window  function  and  a  DFS, 


"  "  t-0 


(6.3) 


where 


Wi£,  f  )  =  ^  5mc{7t/(l-|el)r}  gMN-l-2e)fr  (5.4) 

N  sminJT) 

For  values  of  f<  1/NT,  Equation  (6.4)  is  approximately  equal  to 


iy(e,0)  =(l-lel)r 


(6.5) 


Therefore,  the  expected  value  of  the  cross-ambiguity  function  for  f«0  is  given  by 
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Equation  (6.6)  is  clearly  of  the  form  of  a  discrete  Fourier  series  with  a  triangular 
windowing  fxinction  W(0,z). 


6.2  Properties  of  the  Discrete  Fourier  Series 
For  the  purposes  of  this  study,  the  discrete  Fourier  series  representation  of  the 
function /(e)  will  be  defined  as 


K-l 


Ae)  =  E 


(6.7) 


where  -112  <  e  <  112.  Using  the  orthogonality  property  of  the  complex  exponential, 
the  weights  and  phases  ({)*  are  determined  by  the  following 


1/2 

b^,  e  /(e)  de  . 

The  DFS  representation  of  a  signal  will  now  be  demonstrated  with  an  example. 
For  this  example,  define 


/(e)  =  .56(e)  -  .56(e-.25)  (6.9) 

where  6(e)  is  the  Dirac  delta  function.  Assume  that  the  function  is  to  be  approximated 
with  8  complex  exponentials  {K=8).  Now,  the  weights  b^  and  the  phases  (ji*  can  be 
determined  from  Equation  (6.8)  and  are  given  by 
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(6.11) 


Table  6.1  shows  the  values  of  f(e)  and  fo^sfe)  for  8  discrete  values  of  e.  Notice  that 
the  actual  signal  is  approximated  exactly  for  the  8  match  points  shown.  A  plot  of  the 
series  evaluated  at  128  values  of  e  within  the  defined  interval  is  shown  in  Figure  6.1. 
As  can  be  seen  from  Figure  6.1,  the  DFS  functions  approximate  the  original  function 
exacdy  at  8  match  points.  The  DFS  does  not  match  the  signal  over  the  entire 
interval.  The  reason  that  this  occurs  is  that  only  8  complex  exponentials  are  being 
used  to  represent  a  function  that  is  actually  represented  by  infinitely  many  sinusoids 
because  of  its  infinite  bandwidth.  In  general,  the  DFS  wiU  match  the  original  function 
at  K  points  equally  spaced  between  -112  <  e  <  112  where  K  is  the  number  of 
frequencies  summed  together  to  form  the  estimate.  A  bandlimited  function  could, 
however,  be  represented  exactly. 


6.3.  Matched  FSK/PSK  Signal  Design  Example 
In  this  Section,  an  FSK/PSK  signal  will  be  designed  that  will  be  matched  to 
function  of  Equation  (6.9).  The  expected  value  of  the  cross-ambiguity  function  is 
shown  in  Figure  6.2.  Two  things  should  be  noticed  from  the  Figure.  First,  the  DFS 
has  been  scaled.  Second,  the  mainlobe  of  the  expected  value  of  the  FSK/PSK  cyclic 
cross-correlation  function  has  been  windowed  by  the  triangular  windowing  function. 
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Figure  6.2.  Scaled  DFS  of  Equation  (6.9)  and  the  Expected  Value  of  the  Mainlobe 
of  the  Matched  FSK/PSK  Cross-Ambiguity  Function. 

The  scaling  results  from  the  fact  that  the  Fourier  weights  must  sum  to  imity  in  order  to 
satisfy  the  constraints  of  a  probability  density  function.  Therefore,  an  FSK/PSK  signal 
that  has  an  arbitrary  cross- ambiguity  fimction  mainlobe  can  be  constructed  using  the 
DFS  theory  derived  in  the  previous  section.  Two  sources  of  error  between  the  desired 
result  and  the  realized  result  are  the  effect  of  the  triangular  windowing  function  and 
the  variance  of  the  realized  main  lobe.  The  windowing  is  a  deterministic  functional 
error  and  the  variance  is  a  random  noise-like  error. 

Figures  6.3.a  and  6.3.b  show  a  realization  of  the  mainlobe  of  an  N=64  pulse 
K=8  frequency  Matched  FSK/PSK  cross-ambiguity  function  designed  to  be  matched  to 
the  fimction  of  Equation  (6.9).  Figure  6.3. a  shows  the  real  component  of  the  cyclic 
cross-correlation  function  and  Figure  6.3. b  shows  the  imaginary  component.  Three 
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N=64  /  K=8  Matched  FSK/PSK  Cyclic  Cross-Correlation 


standard  DoviationB  From  Mean  Vdlm  Shown 


Figure  6.3. a.  Real  Component  of  the  Realized  Matched  FSK/PSK  Cross-Ambigmty 
Function  and  Standard  Deviations  for  N=64  and  K=8. 


N=64  /  K=8  Matched  FSK/PSK  Cyclic  Cross-Correlation 


standard  Dmot'ionB  From  Moon  Value  Shown 


Figure  6.3  .b.  Imaginary  Component  of  the  Realized  Matched  FSK/PSK  Cross- 
Ambiguity  Function  and  Standard  Deviations  for  N=64  and  K=8. 
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curves  are  shown  in  each  figure:  the  realized  function,  the  expected  function  plus  the 
standard  deviation  function,  and  the  expected  function  minus  the  standard  deviation 
function.  As  can  be  seen  from  Figures  6.3.a-b,  the  realized  functions  stay  within  the 
previously  derived  standard  deviation  function  where  the  standard  deviation  is  defined 
as  the  square  root  of  the  variance.  Therefore,  Figures  6.3. a  and  6.3.b  verify  the 
expected  value  Equation  (4.5)  and  the  variance  Equation  (4.7). 

The  magnitude  squared  of  the  error  between  the  expected  Matched  FSK/PSK 
cross-ambiguity  function  and  the  actual  realized  Matched  FSK/PSK  cross-ambiguity 
function  of  Figures  6.3. a-b  is  shown  in  Figure  6.4.  As  can  be  seen  from  Figure  6.4, 
the  realized  error  (or  variance)  loosely  follows  the  mainlobe  variance  derived  in  the 
previous  chapter.  Figure  6.5  shows  the  average  realized  variance  of  300  random 
Matched  FSK/PSK  cross-ambiguity  functions  each  for  A=64,  128,  256,  and  512 
pulses.  As  predicted  by  Equation  (4.7),  the  mainlobe  variance  decreases  as  1/N. 
Additionally,  Equation  (4.7)  is  verified  since  the  shape  of  the  variance  curve  is  closely 
tracked  by  the  realized  variance  values.  It  should  be  reiterated  that  Figure  6.4  shows  a 
single  realization  of  a  random  process  (the  error  between  the  realized  mainlobe  and  the 
expected  mainlobe)  while  Figure  6.5  shows  an  average  of  300  realizations.  The  curve 
of  Figure  6.4  is  not  expected  to  follow  exactly  with  the  derived  mean  squared  error 
since  it  is  only  a  single  realization. 

6.4.  Conclusions 

The  mean  value  of  the  cross-ambiguity  function  mainlobe  range  response  is  of 
the  form  of  a  discrete  Fourier  series.  The  weights  of  the  discrete  Fourier  series  are 
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Cyclic  Cross-Correlation  Mainlobe  Variance 


Figure  6.4.  Squared  Error  of  the  Realized  Matched  FSK/PSK  Cross-Ambiguity 
Function  Along  with  the  Expected  Error  Function. 


Figure  6.5.  Average  Realized  Squared  Error  Functions  and  Expected  Error 
Functions  for  N=64,  128,  256,  and  512  Pulses. 
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defined  by  the  probability  density  function  of  the  frequency  selection  process,  and  the 
phases  are  given  by  the  sub-pulse  phases.  Using  known  properties  of  the  discrete 
Fourier  series,  FSK/PSK  signals  can  be  designed  to  approximate  any  arbitrary 
mainlobe  cross-ambiguity  function  range  response.  This  makes  it  possible  to  design 
signals  that  excite  a  maximum  response  from  any  arbitrary  radar  target  based  on 
knowledge  of  the  target  range  profile  h(t).  Additionally,  computer  simulations  verify 
the  results  of  Chapter  IH;  the  variance  of  the  mainlobe  response  decreases  with 
increasing  sequence  length.  The  utility  of  this  signaling  scheme  wiU  be  demonstrated 


in  the  next  chapter. 


CHAPTER  Vn 


COMPARATIVE  ANALYSIS  OF  MATCHED  FSK/PSK 

In  this  chapter,  matched  FSK/PSK  signaling  is  compared  with  two  radar 
signaling  schemes:  a  baseline  system  that  utilizes  a  uniform  spectrum  signal  with  a 
receiver  matched  to  the  transmitted  signal  and  a  matched  filter  system  that  utilizes  a 
post  processing  filter  matched  to  the  target  to  filter  the  output  of  the  baseline  system. 
The  baseline  system  is  referred  to  as  a  current  technology  system  because  its  signal 
processor  is  matched  only  to  the  transmit  signal.  The  matched  filter  system  is  a  future 
technology  system  since  it  makes  use  of  target-specific  scattering  information  (the 
range  profile).  The  analysis  begins  with  an  analysis  of  the  radar- target  interaction 
(effective  radar  cross  section)  of  the  baseline  system,  a  white  noise  contamination 
analysis  of  the  baseline  system,  and  a  clutter  contamination  analysis  of  the  baseline 
system.  Next,  the  effective  target  radar  cross  section,  white  noise  contamination,  and 
ground  clutter  contamination  of  the  matched  filter  radar  is  analyzed  in  a  similar 
fashion.  The  matched  FSK/PSK  radar  system  is  analyzed  for  the  same  three 
conditions  above  to  develop  equations  for  comparison  with  the  other  two  methods.  A 
general  mathematical  comparison  follows  the  analysis  of  the  three  radar  technologies 
which  are  quantified  by  performance  ratios,  or  performance  improvements.  Finally, 
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the  performance  improvements  are  characterized  for  six  example  target  spectra 
predicted  by  a  radar  cross  section  prediction  code. 


7.1.  Analysis  of  the  Baseline  System 

The  baseline  system  transmits  an  FSK/PSK  signal  with  a  frequency  probability 
density  function  defined  by 


(7.1) 


where  K  is  the  number  of  frequencies  used  by  the  radar  signal.  Therefore,  according 
to  Equation  (6.6),  the  ambiguity  function  of  the  radar  is  approximately  given  by 


Xje  .  0)  =  ^  E 

^  t-o 

Equation  (7.2)  has  an  impulse-like  ambiguity  range  delay  function.  It  is  assumed  that 
the  target  has  a  range  profile  defined  by 


*(6) -i:  i/^ 

where  ct*  is  the  radar  cross  section  of  the  target  at  the  A:""  frequency,  and  the  target 
imparts  a  phase  shift  of  on  the  A"*  frequency  component  of  the  signal.  From  the 
correlation  property  of  Equation  (3.49),  the  maximum  possible  output  of  the  radar  for 
a  range  delay  value  of  e=0  is  given  by 
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(7.4) 


It  should  be  noted  that  the  squared  value  of  the  maximum  output  voltage  is  equivalent 
to  the  effective  radar  cross  section  a  of  the  target  used  in  Equation  (2.1). 

The  variance  of  an  additive  white  noise  contaminant  is  given  by  Equation 
(3.36)  to  be 


N 

R  (0)  =  — 


(7.5) 


where  is  the  power  spectral  density  of  the  channel  noise  and  is  the  duration  of 
the  signal.  From  Equation  (4.12),  the  ISL  of  the  signal  with  the  frequency  distribution 
defmed  by  Equation  (7.1)  is  given  approximately  by 


K-l 


ISL  -E  p'-j 


kM) 


(7.6) 


Therefore,  the  output  of  the  baseline  radar  has  been  defined  for  the  effective  radar 
cross  section  (7.4),  additive  noise  (7.5),  and  clutter  contamination  (7.6).  The  clutter 
contamination,  for  the  purposes  of  this  comparative  study,  wUl  be  defined  as  the  ratio 
of  the  effective  target  radar  cross  section  to  the  integrated  sidelobe  level.  Or,  as  the 
ratio  of  the  square  of  Equation  (7.4)  to  Equation  (7.6) 
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1.2.  Matched  Filtering 

A  matched  filter  implementation  is  now  investigated.  A  matched  filter 
implementation  involves  a  radar  that  utilizes  the  baseline  signaling  scheme  with  the 
frequency  probability  density  function  defined  by  Equation  (7.1).  As  in  the  previous 
section,  the  output  of  the  radar  is  given  by  Equation  (7.4).  The  matched  filter  radar 
system  differs  from  the  baseline  system  in  that  it  filters  the  range  profile  (measured  by 
the  baseline  radar)  with  a  filter  impulse  response  defined  by 


AT-l  - 
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(7.8) 
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The  denominator  term  of  Equation  (7.8)  is  a  normalization  factor  to  give  the  filter  a 
unity  noise  bandwidth  [1].  From  linear  system  theory,  it  is  known  that  the  output  of 
the  matched  filter  for  e=0  is  equal  to  the  convolution  of  the  filter  impulse  response  of 
Equation  (7.8)  with  the  cross-ambiguity  function  of  Equation  (7.4) 
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(7.9) 


As  above  with  Equation  (7.4),  the  squared  value  of  the  output  voltage  is  equal  to  the 
effective  target  radar  cross  section  of  the  target  a.  It  should  be  noticed  that  Equation 
(7.9)  does  not  contain  any  phase  terms.  This  results  from  the  fact  that  the  matched 
filter  is  phase  matched  to  the  target,  and  therefore,  realigns  all  of  the  phase  variations 
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in  the  target  spectrum.  In  other  words,  all  contributions  to  the  target  spectrum  are 
added  in  phase  with  one  another  (added  constructively). 

For  the  baseline  radar,  the  autocorrelation  function  of  the  noise  process  at  the 
input  of  the  matched  filter  is  given  by  Equation  (3.36)  to  be 


N  I  .. 
R  (e)  =  —  T  —  eJ' 

n/1  '  '  rn  1/- 
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(7.10) 


From  stochastic  processes  theory  [1],  the  autocorrelation  function  of  the  noise  process 
output  from  the  matched  filter  defined  by  Equation  (7.10)  is  given  by 
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The  mean  squared  noise  output  voltage  (the  noise  power)  output  from  the  matched 
filter  is  given  by 


0 )  = 
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K  r 


(7.12) 


Therefore,  the  noise  variance  of  the  baseline  radar  (7.5)  is  K  times  that  of  the  matched 
filtered  case. 

To  determine  the  effect  of  the  matched  filter  on  the  clutter  return,  it  is  assumed 
that  the  clutter  can  be  modeled  as  additive  white  noise.  In  other  words,  the 


autocorrelation  function  of  the  clutter  return  can  be  modeled  according  to  (7.11)  as 
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where  the  ISL  is  given  in  Equation  (7.6).  Following  the  white  noise  analysis 
procedure  from  above,  the  actual  effect  of  the  clutter  variance  is  equal  to 

R  (0)^—-—  .  (7.14) 

“  K 

Therefore,  the  signal  to  clutter  ratio  for  the  matched  filter  case  is 
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7.3.  Matched  FSK/PSK  Radar 

The  matched  FSK/PSK  radar  is  said  to  be  matched  to  a  target  if  the  probability 
density  function  of  the  frequency  selection  process  has  a  shape  that  is  proportional  to 
the  magnitude  of  the  target’s  voltage  reflection  coefficients  at  each  of  the  frequencies. 
In  equation  form,  the  probability  density  function  is  given  by 


iW) 


(7.16) 


The  numerator  of  Equation  (7.16)  describes  the  shape  of  the  probability  density 
function  and  the  denominator  is  a  normalizing  factor  required  of  a  valid  probability 
density  function.  The  cross-ambiguity  function  for  the  matched  FSK/PSK  radar 
waveform,  as  defined  by  Equation  (6.6),  is  therefore  equal  to 


82 


X,(e,0) 

where  the  Pi^’s  are  as  defined  above  in  Equation  (7.16)  and  the  phase  terms  (j)^^  are  also 
given  in  Equation  (7.16).  It  is  of  interest  to  determine  the  value  of  the  output  signal 
for  zero  delay.  In  other  words,  when  the  cross-ambiguity  function  is  perfectly  aligned 
with  the  target  range  profile  of  Equation  (7.3).  Therefore,  according  to  Equation 
(3.49),  the  peak  response  for  e=0  equals 
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As  with  the  previous  two  radar  examples.  Equation  (7.18)  is  the  square  root  of  the 
effective  target  radar  cross  section  of  the  target  a. 

It  is  now  of  interest  to  investigate  white  noise  contamination  of  a  matched 
FSK/PSK  radar.  According  to  Equation  (3.36),  for  the  case  of  the  matched  FSK/PSK 
radar,  noise  autocorrelation  function  equals 
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Tj  k^ 


where  the  p^’s  are  as  given  above  in  Equation  (7.16),  Ng  is  the  channel  noise  power 


spectral  density,  and  Tj  is  the  waveform  duration.  For  the  case  of  e=0,  the  mean 


noise  output  power  is 
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(0)  =  ^ 


(7.20) 


since  the  probabilities  sum  to  unity.  Notice  that  this  is  the  same  as  for  the  baseline 
case  above  (7.5). 

From  Equation  (4.12),  the  ISL  of  a  matched  FSK/PSK  radar  is  given 
approximately  by 
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Therefore,  the  signal  to  clutter  ratio  (as  defined  by  the  ratio  of  effective  radar  cross 
section  of  the  target  to  integrated  sidelobe  level)  is  given  by 


5,  _  d(  0  ,  0  )  _  ^^  ^  (7.22) 

c  ISL  ri  '  ' 

Notice  that  the  signal  to  clutter  ratios  of  the  matched  FSK/PSK  radar  (7.22)  and  the 
matched  filter  radar  (7.15)  are  identical. 


7.4.  Comparison  of  Results 

In  this  section,  the  results  obtained  in  the  previous  three  sections  are  compared. 
The  signaling  schemes  are  compared  in  three  categories:  effective  target  radar  cross 
section  (radar  cross  section  enhancement),  additive  noise  contamination  level,  and 


signal  to  clutter  ratio. 
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A  gain  can  be  defined  that  is  the  ratio  of  the  matched  FSE7PSK  effective  target 
radar  cross  section  obtained  from  Equation  (7.18)  to  the  effective  target  radar  cross 
section  of  the  baseline  radar  from  Equation  (7.4).  This  ratio  describes  the  target 
enhancement  obtained  by  utilizing  matched  FSK/PSK  signaling.  It  should  be  noted 
that  this  is  a  power  gain  and  therefore  translates  into  a  signal  to  noise  ratio 
improvement.  The  target  enhancement  of  matched  FSK/PSK  is  therefore  given  by 
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It  can  be  seen  that  Equation  (7.23)  is  always  greater  than  or  equal  to 


(7.24) 


Notice  that  the  enhancement  is  target  dependent  because  the  matched  FSK/PSK  radar 
system  utilizes  target-specific  information  while  the  baseline  system  does  not. 

A  gain  can  also  be  defined  that  is  the  ratio  of  the  matched  FSK/PSK  effective 
target  radar  cross  section  to  the  matched  filter  target  radar  cross  section.  This  is  given 
by  the  squared  ratio  of  Equation  (7.18)  to  Equation  (7.9)  and  is  equal  to 
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There  is  considered  to  be  a  performance  improvement  due  to  matched  FSK/PSK 
signaling  over  a  matched  filter  radar  if 


(7.26) 


Using  the  triangle  inequality,  the  left-hand  side  of  Equation  (7.26)  represents  the 
hypotenuse  while  the  right-hand  side  represents  the  average  side  length.  Therefore, 
Equation  (7.26)  is  satisfied  since  the  hypotenuse  is  always  greater  than  the  average 
side  length.  Therefore,  the  effective  radar  cross  section  is  always  greater  for  the 
Matched  FSK/PSK  radar  than  for  the  matched  filtered  radar. 

Target  enhancements  for  two  example  target  radar  cross  section  distributions 
can  now  be  calculated:  a  uniform  spectral  distribution  and  an  impulsive  spectral 
distribution.  For  the  uniform  distribution,  =  c  and  (])*  =  0  for  all  k  from  0  to  ^-1 
where  K  is  the  number  of  frequencies  in  the  signal.  For  this  specific  case,  the 
matched  FSK/PSK  target  enhancement  to  baseline  radar  target  enhancement  is  given 


by 
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£  =  1  (7.27) 

D 

The  enhancement  of  the  two  radar  targets  is  identical  since  the  target  spectrum  is 
identical  to  the  baseline  radar  transmit  signal  spectrum.  In  other  word,  the  baseline 
radar  is  matched  to  the  uniform  spectrum  target.  Likewise,  the  enhancement  of 
matched  FSK/PSK  radar  over  matched  filter  radar  is 

E  =K  .  (7.28) 

m 

The  matched  filter  radar  suffers  a  loss  in  target  enhancement  for  the  impulsive  target 
even  though  the  transmit  signal  is  matched  to  the  target.  This  is  a  result  of  the  unity 
noise  bandwidth  of  the  filter.  In  other  words,  the  filter  tends  to  attenuate  the  target 
return,  but  as  will  be  seen  later,  the  noise  is  attenuated  by  the  same  amount.  For  the 
impulsive  target  distribution  the  target  spectrum  equals  o„  =  c  for  some  integer  m  and 
is  zero  otherwise.  For  this  impulsive  spectral  distribution,  the  matched  FSK/PSK 
target  enhancement  over  the  baseline  radar  is  equal  to 

E=K^  (7.29) 

O 

The  target  enhancement  for  this  spectral  distribution  is  quite  large.  The  reason  for  this 
improvement  is  that  the  baseline  signal  only  has  one  frequency  component  reflected 
whereas  the  matched  FSK/PSK  signal  only  transmits  in  the  frequency  band  reflected 
by  the  target.  The  target  enhancement  of  matched  FSK/PSK  over  matched  filtering  is 


given  by 
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E  .  (7-30) 

m 

Therefore,  the  target  enhancement  is  dependent  upon  the  distribution  of  the  target 
radar  cross  section  in  frequency.  The  target  enhancement  is  greater  for  the  impulsive 
spectral  distribution  (uniform  in  range)  than  for  the  uniform  spectral  (impulsive  in 
range)  radar  cross  section  distribution.  This  is  a  result  of  the  effective  utilization  of 
the  transmitter  power  for  the  matched  FSK/PSK  radar. 

From  the  noise  analysis  of  the  previous  sections,  the  improvement  in  signal  to 
noise  ratio  of  the  matched  FSK/PSK  radar  SNR/^  over  the  baseline  radar  SNR^  is 
simply  given  by  Equation  (7.23)  since  the  noise  contamination  levels  are  the  same  for 
each  of  the  two  systems.  From  Equations  (7.12),  (7.20),  and  (7.25),  the  improvement 
in  signal  to  noise  ratio  of  the  matched  FSK/PSK  radar  system  SNR/^  over  that  of  the 
matched  filter  radar  SNR„  is  given  by 
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(7.31) 


For  the  two  example  distributions  from  above,  the  signal  to  noise  ratio  improvements 
are 


SNR^ 

SNR, 


(7.32) 


and 
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for  the  uniform  spectral  distribution,  and 


(7.33) 


(7.34) 

SNR^  “ 
and 

^  =  El  =  K  (7.35) 

SNR  K  K 

m 

for  the  impulsive  spectral  distribution.  This  is  intuitively  pleasing  since  the  uniform 
spectral  distribution  matches  the  signal  transmitted  of  the  baseline  implementation  (in 
other  words,  the  target-specific  information  is  the  same  for  all  targets).  Therefore,  no 
performance  improvement  should  be  expected.  The  performance  improvement 
associated  with  the  impulsive  spectral  distribution  is  also  logical  since  the  radar  would 
be  transmitting  only  in  the  spectral  region  in  which  the  target  reflects. 

The  signal  to  clutter  performance  of  the  matched  FSK/PSK  system  (7.22)  is 
identical  to  that  of  the  matched  filter  system  (7.15).  It  should  be  remembered  that  the 
matched  filtering  implementation  has  an  additional  layer  of  signal  processing  over  the 
matched  FSK/PSK  radar  and  hence  is  more  complex  in  the  receiver  whereas  matched 
FSK/PSK  is  more  complex  in  the  transmitter.  The  matched  FSK/PSK  radar  system 
requires  that  the  signal  be  transmitted  and  processed  through  the  correlation  processor. 
The  matched  filter  implementation  requires  that  the  output  of  the  correlator  be 
matched  filtered  to  enhance  the  target. 
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The  ratio  of  the  matched  FSK/PSK  signal  to  clutter  ratio  to  the  baseline  signal 
to  clutter  ratio  is  equal  to 
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(7.36) 


From  Equations  (7.27)  and  (7.29),  it  is  known  that  Equation  (7.36)  varies  between  one 
for  the  uniform  spectral  distribution  and  K  for  the  impulsive  spectral  distribution. 
Therefore,  either  matched  FSK/PSK  signaling  or  matched  filtering  provide  a  gain  over 
the  baseline  radar  in  all  environments  due  to  the  use  of  target-specific  information. 


7.5.  Example  Performance  Gains  for  Predicted  Target  Signatures 
The  calculations  in  the  previous  examples  are  for  some  extreme  theoretical 
distributions.  It  is  of  interest  to  calculate  the  performance  improvements  for  some 
predicted  signatures  of  actual  aircraft  models.  Table  7.1  shows  the  performance 
improvements  in  dB  for  six  different  aircraft  targets.  The  target  spectra  (a^’s  and  (jilt’s) 
are  predicted  for  the  models  using  Xpatchl  [2]  which  is  an  optical  region 
electromagnetic  signature  prediction  code.  The  signatures  are  taken  over  an  azimuthal 
sweep  from  5°  to  15°  off  nose-on  in  10  equally  spaced  steps  and  in  elevation  from  -5° 
to  5°  above  the  waterline  in  10  equally  spaced  steps.  There  are  a  total  of  100 
signatures  for  each  target.  Each  signature  is  predicted  within  a  bandwidth  from  2.0  to 
2.3  GHz  (1/2  meter  range  resolution),  and  there  are  64  frequencies  in  each  signature. 
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In  Table  7.1,  the  first  column  represents  the  mean  value  of  the  improvement  in 
signal  to  additive  noise  power  of  a  matched  FSK/PSK  radar  over  a  baseline  radar  in 
dB.  This  is  calculated  by  the  following  equation 


=  10  log,,  (£{£,}) 

where  E{}  is  the  ensemble  average  over  aU  azimuth  and  elevation  values,  and  is 
calculated  from  Equation  (7.23)  for  each  azimuth  and  elevation  value. 

Column  two  represents  similar  values  of  improvement  in  signal  to  noise  ratio 
for  the  matched  FSK/PSK  radar  over  the  matched  filter  radar.  The  equation  evaluated 
in  column  two  is  given  by 


=  10  logj„ 

where  E^IK  is  calculated  using  Equation  (7.31). 

Column  three  represents  the  improvement  in  the  signal  to  clutter  ratio  by  using 
matched  FSK/PSK  or  matched  filtering  over  the  baseline  radar  system.  Column  three 
is  calculated  using  the  following  equation; 
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where  the  ratio  in  the  argument  of  the  logarithm  is  calculated  using  Equation  (7.36). 
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Table  7.1.  Realized  Performance  Improvements  and  Losses  for  Matched  FSK/PSK 
and  Matched  Filtering. 


Target  # 

E, 

E. 

S/C 

ElPI 

1 

6.2 

0.9 

5.3 

-3.1 

2 

4.2 

0.8 

3.3 

-2.7 

3 

3.0 

0.5 

2.4 

-2.1 

4 

5.6 

0.9 

4.7 

-3.2 

5 

3.6 

0.8 

2.8 

-2.7 

6 

3.6 

0.7 

2.8 

-2.5 

Units 

dB 

dB 

dB 

dB 

Column  four  of  Table  7.1  shows  the  expected  LPI  performance  losses  due  to 
the  reduction  in  effective  bandwidth  that  results  for  utilizing  matched  FSK/PSK 
signaling.  This  loss  is  calculated  using  Equation  (5.19)  and  is  given  by 

L^,  ^  10  loglO  {  E  {  L^,  ]  )  .  (5.40) 

The  data  in  Table  7.1  show  that  the  performance  improvement  achieved  by 
making  use  of  target-specific  information  is  significant.  For  the  matched  FSK/PSK 
signal  to  noise  ratio  data  (column  one),  the  mean  improvement  over  the  baseline  radar 
is  in  a  range  of  between  3  dB  and  6  dB.  In  other  words,  the  average  signal  to  noise 
ratio  is  doubled  or  even  quadrupled  depending  on  the  target.  For  clutter  limited 
environments,  the  mean  signal  to  clutter  ratios  (column  three)  are  improved  by  a  factor 
of  between  2.4  dB  and  5.3  dB.  This  improvement  is  not  quite  as  large  as  those  seen 
for  the  additive  noise  environments,  but  they  are  still  significant  improvements.  The 
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data  in  column  two  show  that  the  mean  improvement  realized  by  utilizing  matched 
FSK/PSK  instead  of  matched  filtering  is  minimal.  They  range  from  a  0.5  dB 
improvement  to  0.9  dB.  Column  four  shows  the  LPI  performance  loss  due  to  the 
reduction  in  effective  bandwidth  due  to  matched  FSK/PSK.  It  appears  that  typical 
losses  are  in  a  range  of  2  dB  to  3  dB.  Therefore,  the  overall  LPI  performance 
improvement  achieved  by  utilizing  matched  FSK/PSK  radar  over  the  baseline  radar  is 
in  a  range  of  between  1  dB  and  3  dB  on  the  average. 

Figures  7.2  and  7.2  give  some  insight  into  the  enhancement  and  loss  values 
encountered  in  Table  7.1.  Figure  7.1  shows  a  typical  range  profile  and  spectrum  for 
Target  1  while  Figure  7.2  shows  a  typical  range  profile  and  spectrum  for  Target  3. 
From  Table  7.1,  Target  1  is  seen  to  have  the  highest  enhancement  and  signal  to  clutter 
improvement  of  all  targets.  Conversely,  Target  3  has  the  worst  enhancement  and 
signal  to  clutter  improvement.  The  reason  for  these  results  can  be  seen  by  comparing 
Figure  7.1  with  Figure  7.2.  The  range  profile  of  Target  1  is  more  distributed  in  range 
(spectrum  is  more  concentrated  in  frequency)  than  the  profile  of  Target  3.  This  means 
that  greater  target  enhancement  should  occur  for  Target  1  since  its  spectrum  is  less 
uniform  in  frequency  (more  distributed  in  range).  Additionally,  it  should  be  noted  that 
the  spectra  of  both  targets  are  broadbanded  and  uniform.  This  is  the  reason  the 
enhancements  and  signal  to  clutter  ratio  improvements  are  closer  to  unity  than  K. 

This  spectral  uniformity  is  a  characteristic  of  optical  region  target  scattering  (impulsive 
in  range,  broadband  in  frequency)  and  therefore,  large  gains  (>  10  dB)  should  not  be 
expected  for  optical  region  targets. 
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1.6.  Conclusions 

Matched  FSK/PSK  radar  and  the  matched  filter  radar  both  yield  performance 
that  is  superior  to  the  baseline  radar  based  on  an  analysis  of  the  signal  to  noise  ratio 
and  the  signal  to  clutter  ratio.  The  performance  increase  is  dependent  upon  the  target 
spectral  content  since  the  matched  FSK/PSK  radar  and  the  matched  filter  radar  both 
make  use  of  target-specific  information.  The  theoretical  matched  FSK/PSK  signal  to 
noise  ratio  is  between  1  and  times  that  of  the  baseline  radar.  Additionally,  the 
matched  FSK/PSK  signal  to  clutter  ratio  is  theoretically  between  1  and  K  times  that  of 
the  baseline  radar  depending  on  the  spectrum  of  the  target.  Therefore,  a  significant 
performance  increase  can  potentially  be  obtained  by  utilizing  target-specific 
information  in  the  design  of  a  radar  signal. 

Based  on  the  calculations  in  Section  7.5,  expected  values  of  signal  to  noise 
ratio  improvement  for  actual  optical  region  targets  are  typically  between  3  dB  and  6 
dB.  Additionally,  the  use  of  target-specific  information  results  in  typical  signal  to 
clutter  improvements  of  between  2  dB  and  5  dB.  Even  though  these  improvements 
are  significant,  they  are  much  less  than  the  theoretically  possible  improvements.  This 
is  a  result  of  the  basically  uniform  spectral  signature  of  optical  region  radar  targets. 
The  spectral  signatures  have  a  uniform  broadband  shape  because  of  the  impulse-like 
temporal  characteristic  of  optical  region  radar  targets. 

In  an  additive  white  noise  environment,  matched  FSK/PSK  radar  outperforms 
both  the  baseline  system  and  the  matched  filter  system.  This  performance  increase 
results  from  an  efficient  use  of  transmitter  power.  The  matched  FSK/PSK  radar 
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transmits  power  only  in  regions  that  scatter  efficiently.  In  other  words,  a  radar  that 
transmits  an  impulse-like  signal  wastes  power  in  spectral  regions  that  do  not  reflect 
transmitter  power.  The  overall  performance  gain  is  dependent  upon  the  target 
reflection  characteristics,  and  the  gain  is  bounded  between  1  and  K  where  K  is  the 
number  of  frequencies. 

The  results  of  calculations  based  on  predicted  target  signatures  show  that  the 
realized  matched  FSK/PSK  signal  to  noise  ratio  improvement  relative  to  the  matched 
filter  radar  is  typically  less  than  1  dB.  Therefore,  it  seems  that  matched  FSK/PSK 
performance  is  not  signiflcantly  better  than  matched  filtering  for  a  actual  optical  region 
target.  Additionally,  with  the  matched  filter  implementation,  there  is  no  LPI 
performance  degradation  due  to  a  transmit  bandwidth  reduction. 

In  a  ground  clutter  limited  environment,  matched  FSK/PSK  radar  performance 
is  identical  to  matched  filtering  radar.  This  results  from  the  fact  that  the  ISL  of  the 
matched  FSK/PSK  radar  signal  increases  identically  with  the  target  enhancement  of 
matched  FSK/PSK  over  matched  filtering.  It  should  be  noted,  that  matched  FSK/PSK 
requires  one  less  layer  of  receiver  signal  processing  than  matched  filtering  because  the 
matched  FSK/PSK  signaling  scheme  has  the  additional  layer  of  signal  processing  built 
into  the  waveform.  This  post  processing  could  be  a  limiting  factor  if  a  large  range 
search  (large  number  of  range  bins)  is  required  because  of  the  computational  overhead 
of  performing  the  discrete  convolutions.  According  to  [3],  an  point  convolution  (N^ 
is  the  number  of  output  range  bins)  requires  approximately 
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3N  logj  iV,  (7-43) 

complex  floating  point  operations  when  an  efficient  FFT-based  convolution  routine  is 
used.  For  example,  a  1024  point  output  result  requires  approximately  31,000  complex 
floating  point  operations.  This  translates  into  775  psec  of  processing  time  for  a  40 
Mflop  processor.  It  should  be  remembered  that  the  775  psec  is  required  for  the 
convolution  alone.  This  does  not  include  any  other  computer  tasking  that  might  be 
required  such  as  Doppler  processing. 

With  regard  to  the  overall  LPI  performance,  the  matched  FSK/PSK  radar  trades 
2  dB  to  3  dB  in  LPI  performance  in  order  to  obtain  an  additional  3  dB  to  6  dB  in 
signal  to  noise  ratio  performance.  This  is  a  positive  overall  LPI  improvement  of 
between  1  dB  and  2  dB  over  the  baseline  radar.  The  reason  that  this  loss  occurs  is 
that  the  matched  FSK/PSK  radar  must  reduce  its  equivalent  transmitter  bandwidth  to 
transmit  only  in  highly  reflective  regions.  Therefore,  the  equivalent  time-bandwidth 
product  is  reduced. 

It  should  be  reiterated  that  both  the  matched  FSK/PSK  radar  and  the  matched 
filter  radar  perform  better  than  the  baseline  radar  in  both  white  noise  environments  and 
ground  clutter  limited  environments.  This  performance  improvement  is  a  result  of 
target  enhancement.  This  target  enhancement  requires  knowledge  about  the  target’s 
scattering  characteristics.  Since  this  knowledge  is  typically  not  known  a  priori,  a 
technique  for  determining  the  target  spectral  characteristics  must  be  developed.  This 
is  the  focus  of  the  remainder  of  this  study. 
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CHAPTER  Vm 


GRADIENT  DESCENT  METHOD  FOR  DESIGNING 
OPTIMUM  RADAR  SIGNALS 

In  previous  chapters,  the  pulse  compression  radar  processing  technique 
designated  matched  FSK/PSK  was  analyzed.  This  processing  technique  makes  it 
possible  to  design  signals  with  compressed  main  lobes  of  arbitrary  pre-designed 
shapes.  Since  the  response  of  a  radar  to  a  given  target  is  given  by  the  correlation  of 
the  compressed  main  lobe  with  the  target  range  profile  (for  a  high  range  resolution 
radar  signal),  it  is  possible  to  design  a  signal  that  will  excite  an  optimum  response 
from  a  target.  A  problem  arises  with  this  method  because  the  design  of  a  signal  to 
excite  the  optimum  response  requires  a  priori  knowledge  of  the  target.  Typically,  the 
range  profile  will  not  be  known  a  priori  and  therefore,  it  is  desirable  to  have  a  method 
for  iteratively  designing  the  optimum  signal  based  on  previous  radar  returns  in  real 
time.  This  chapter  focuses  on  a  method  that  is  derived  based  on  the  technique  of 
gradient  descent. 

First,  a  brief  outline  of  the  required  background  information  is  given  so  that  the 
problem  can  be  properly  formulated.  Next,  the  optimum  signal  is  derived  for  a 
randomly  varying  target  range  profile.  An  iterative  solution  to  the  analytic  problem  is 
then  derived  that  is  based  on  gradient  descent  [1].  Algorithm  implementation  within 
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the  framework  of  a  radar  is  discussed  next.  Finally,  the  results  of  computer 
simulations  are  used  to  illustrate  the  effectiveness  of  the  algorithm  on  high  range 
resolution  radar  targets  in  various  environments;  additive  white  noise,  correlated  and 
imcorrelated  random  range  variation,  and  correlated  and  uncorrelated  random  angle 
variation.  Conclusions  about  the  usefulness  of  the  algorithm  as  a  target  learning,  or 
measurement  technique  are  discussed. 

8.1.  Background 

A  block  diagram  of  a  high  range  resolution  radar  system  utilizing  a  cross¬ 
correlation  receiver  is  shown  in  Figure  3.1.  The  response  of  a  high  range  resolution 
radar  range  profile  h(t)  to  a  radar  utilizing  a  cross-correlation  receiver  is  defined  by 


-00 

where  the  cross-ambiguity  function  XxyW  given  by  Equation  (6.6). 

An  optimum  radar,  or  matched  radar  is  defined  as  a  radar  that  utilizes  a  cross¬ 
ambiguity  that  maximizes 


00  oo 

Xj  0  ,  0  )p  =  J  J  x’^(y,0)  h{u)h\v)du  dv 


(8.2) 


It  is  useful  to  cast  the  equation  above  into  a  vector  representation.  The  vector 
representation  of  the  cross-ambiguity  function  is 
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X=[Xj^o)  . 

where  the  subscripts  designate  the  various  sample  points  and  T  designates  the  vector 
transposed.  Likewise,  for  the  target  range  profile  is  given  by 


h  =[/i(Mo)  . hiu^_^)Y 

Equation  (8.2)  becomes 


(8.4) 


k(0,0)p  =  x»/»/i"x 

where  the  superscript  H  designates  the  complex  conjugate  transpose.  Therefore,  it  is 
desired  to  maximize  Equation  (8.5)  with  respect  to  x-  For  generality,  it  is  assumed 
that  the  vector  h  is  a  random  process  described  by  a  correlation  matrix  C  defined  by 


C^E{hh^} 

where  E  represents  the  ensemble  average  operator.  Therefore,  Equation  (8.5)  becomes 


lx(o,o)p  =  x^Cx 

where  C  is  the  target  correlation  matrix. 


(8.7) 


8.2.  Analytic  Approach  to  Optimization 
Equation  (8.7)  must  now  be  maximized.  It  is  of  interest  to  maximize  the 
squared  voltage  value  (as  opposed  to  the  signal-to-noise  ratio)  since  the  white  noise 
contamination  of  a  pulse  compression  waveform  depends  only  on  the  processing 
interval  and  not  the  compressed  main  lobe  shape.  The  typical  approach  to  maximizing 
a  functional  is  to  take  the  derivative  with  respect  to  the  control  function,  and  then 
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solve  for  the  value  of  the  control  function  for  which  the  derivative  is  zero.  It  is  also 
required  that  the  energy  of  the  cross-ambiguity  function  be  constrained  to  some 
constant  value  a.  This  is  a  constrained  extremization  problem,  therefore  a  Lagrange 
multiplier  must  be  used.  Using  the  Lagrange  multiplier,  the  following  equation  must 
be  maximized 

x„ix) 

where  X  is  known  as  the  Lagrange  multiplier.  The  optimum  solution  to  Equation  (8.8) 
is  obtained  by  solving 


9)C„(X) 

^X 


=2Cx-2Xx=0 


(8.9) 


Therefore,  the  optimum  vector  Xopt  is  obtained  from  the  following 


Equation  (8.10)  is  simply  a  statement  of  the  well-known  eigenvalue  problem  where 
Xop,  is  an  eigenvector  of  C  and  X  is  its  associated  eigenvalue.  It  can  be  shown  that  the 
optimum  solution  of  Equation  (8.10)  occurs  when  Xopt  is  the  eigenvector  of  C 
associated  with  the  principal  eigenvalue  A  of  C.  x<,p,  is  also  known  as  the  principal 
component  of  the  correlation  matrix  C.  The  principal  component  can  be  thought  of  as 
pointing  in  the  direction  of  the  maximum  variation  of  the  random  process  h  [2]. 

Equation  (8.10)  suggests  how  to  design  a  signal  based  on  a  target  correlation 
matrix.  This  correlation  matrix  is  typically  not  known  a  priori.  Therefore,  an  iterative 
solution  must  be  obtained  that  will  iteratively  approximate  the  optimum  solution. 
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8.3.  Gradient  Descent  Iterative  Solution 
An  iterative  method  for  maximizing  Equation  (8.8)  is  described  in  this  section. 
First  an  initial  guess  of  Xa  is  assumed.  Also,  it  is  assumed  that  all  subsequent 
estimates  of  the  optimum  solution  will  be  better  than  the  previous  estimate.  In 
equation  form,  the  previous  assumption  yields 

l.A'U)  ^  (  Xo  + 

Expanding  x*r  about  Xo  as  a  Taylor  series  and  retaining  only  the  linear  term  yields  the 
following  linearized  equation 


tr  (  Xo  +  ^X  )  «  tr  (  Xo  )  +  ^X 


9  X, 


(8.12) 


Substituting  Equation  (8.12)  into  Equation  (8.11)  gives  the  following  inequality 


Xxr  (  Xo  )  ^  X„  (  Xo  )  +  Ax  . 

Equation  (8.13)  is  always  satisfied  if  the  increment  Ax  is  given  by 


(8.13) 


Ax  =  P 


^  Ixr 

d  X 


where  p  is  the  learning  rate.  Therefore, 


X  =  X  *  +  P 


dR 
d  X 


(8.14) 


(8.15) 


where  the  subscript  designates  the  iteration  index.  Using  the  partial  derivative  of  the 
functional  from  the  previous  section.  Equation  (8.15)  becomes 
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Xt.i  =X*  +  2y 

Another  crucial  component  of  the  iterative  solution  involves  normalizing  Xk+i 
each  step  with  the  following  energy  normalizing  equation 


X 


k  +  1 


X  t  +  1 


^X  t  + 1  X  *  + 1 


(8.17) 


This  normalization  process  is  required  since  the  optimization  problem  is  a  constrained 
optimization. 

It  is  important  to  note  that  the  algorithm  defined  by  Equations  (8.16)  and  (8.17) 
is  known  as  plain  Hebbian  Learning  [2]  in  the  field  of  Neural  Sciences.  The  learning 
technique  can  be  thought  of  as  a  reinforced  learning  technique  in  that  it  tends  to 
converge  to  the  vector  that  points  in  the  direction  of  maximum  variation  (variance).  It 
will  also  point  in  the  direction  of  the  vector  input  that  is  introduced  to  the  algorithm 
most  often.  The  learning  machine  paradigm  of  the  iterative  design  algorithm  yields 
significant  intuitive  insight  into  the  operation  of  the  algorithm. 


8.4.  Implementation  Issues 

Equation  (8.16)  provides  a  method  for  designing  the  optimum  radar  cross¬ 
ambiguity  function.  The  technique  requires  the  target  correlation  matrix  C  to  be 
measured.  Therefore,  a  technique  for  measuring  C  must  be  determined.  From 
Equation  (8.6), 
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C=E{hh^} 
and 

Cx  =£(**"}  X 

Now,  define 

Cx-»(*”x) 

^  h  c- 

where  c  is  a  scalar  that  is  a  measure  of  the  cross-correlation  between  h  and  %■  The 
radar  can  measure  c  by  transmitting  the  complex  conjugate  of  %  and  conjugating  the 
received  voltage  for  zero  delay,  h  can  be  measured  if  the  radar  transmits  a  signal  with 
an  impulse-like  cross -ambiguity  function.  Therefore,  a  plausible  implementation  of  the 
gradient  descent  algorithm  derived  in  the  previous  section  would  be  the  following: 

1)  Transmit  the  complex  conjugate  of  Xk  and  receive  c, 

2)  Transmit  an  impulse  like  waveform  and  receive  h,  and 

3)  Form  the  next  estimate  of  the  optimum  waveform  from 

Xfc.i=X*-"2pc*A  (8.21) 

The  algorithm  above  requires  that  the  optimum  waveform  be  transmitted  only 
half  of  the  time.  The  other  half  of  the  time,  an  impulse-like  sub-optimum  signal  must 
be  transmitted.  This  is  required  only  during  the  waveform  design  procedure.  In  other 
words,  the  radar  would  operate  in  a  design  mode  until  the  signal  is  designed  (i.e.,  K 


(8.18) 


(8.19) 
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iterations).  After  the  iteration,  the  radar  would  operate  using  only  the  designed 
optimum  waveform.  This  aspect  is  illustrated  in  the  following  computer  simulations. 

8.5.  Computer  Simulation  of  the  Gradient  Descent  Method 
In  this  section,  the  gradient  descent  signal  design  technique  derived  in  the 
previous  sections  is  illustrated  using  computer  simulations.  These  simulations  will 
model  the  effects  of  additive  white  noise,  correlated  and  uncorrelated  random  range 
motion,  and  correlated  and  uncorrelated  random  angular  motion.  A  listing  of  the 
Matlab®  [3]  script  used  in  the  study  can  be  found  in  Appendix  A. 

The  target  range  profile  used  in  the  simulations  is  that  of  an  aircraft-like  target 
generated  by  Xpatchl  [4].  The  range  profiles  consist  of  64  data  points  and  they  are 
valid  within  the  bandwidth  of  2.0  to  2.3  GHz.  Therefore,  they  have  a  .5  meter 
resolution.  The  angular  variation  is  between  5°  and  15°  in  azimuth  at  10°  below  the 
waterline  in  elevation.  These  angular  locations  were  chosen  since  they  are  typical  of  a 
closing  target.  Additionally,  it  was  found  that  the  incident  angle,  at  least  on  a  broad 
scale,  had  little  to  do  with  the  performance  of  the  algorithm. 

The  performance  metric  is  the  squared  error  between  the  estimated  range 
profile  (optimum  signal)  and  the  actual  range  profile 


et  =  iX,  -  {%,  -  h) 


(8.22) 


where 
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X/  h  (8.23) 

is  the  cross-correlation  between  the  estimated  profile  and  the  actual  range  profile.  It 
has  been  assumed  that 

z“x,  -»"<>  =  1 

The  actual  squared  error  value  in  Equation  (8.22)  would  replace  the  magnitude 
operator  with  an  operator  that  would  include  the  real  part  only.  The  magnitude 
operator  is  used  for  this  case  since  the  phase  error  is  not  of  interest,  only  the  shape. 
Also  notice  that  the  iterative  algorithm  maximizes  the  component  of  Equation  (8.23) 
and  therefore  Equation  (8.22)  is  minimized.  Therefore,  the  gradient  descent  algorithm 
minimizes  the  squared  error  between  the  optimum  signal  and  the  actual  target  range 
profile.  This  means  that  the  gradient  descent  technique  is  optimal  for  target 
measurement  as  well  as  target  excitation. 

8. 5. A.  Additive  White  Noise  Effects 

Figure  8.1  shows  the  effect  of  the  learning  rate  on  the  removal  of  additive 
white  Gaussian  noise.  This  plot  shows  the  mean  square  error  versus  the  number  of 
iterations  for  three  values  of  learning  rate:  fi=0.2,  0.1,  and  0.05.  The  average  of  100 
independent  trials  is  shown  for  each  value  of  the  learning  parameter  ji.  The  signal  to 
noise  ratio  of  each  measurement  is  unity  where  the  signal  to  noise  ratio  is  defmed  by 
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Additive  White  Noise  Effects 


Figure  8.1.  Mean  Squared  Error  Versus  the  Number  of  Iterations  for  Three 
Different  Learning  Rates. 


SNR  =  (8.25) 

n"n 

where  n  is  the  vector  of  64  independent  Gaussian  random  variables.  From  the  plot  it 
can  be  seen  that  the  gradient  descent  technique  converges  to  its  minimum  value  after 
1/ji  iterations,  and  the  mean  square  error  converges  to  approximately  a  value  of  ji. 
Figure  8.2  shows  the  effect  of  additive  noise  for  a  fixed  learning  rate.  The  plot  shows 
the  mean  square  error  for  three  different  signal  to  noise  ratios:  SNR  =  6  dB,  0  dB,  and 
-6  dB.  The  learning  rate  is  a  constant  value  of ;/  =  0.1.  An  average  of  100 
independent  trials  is  shown  for  each  curve.  It  can  be  seen  from  the  Figure  that  the 
mean  square  error  of  the  estimate  improves  at  10  dB/decade  with  signal  to  noise  ratio. 
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Additive  White  Noise  Effects 


Figure  8.2.  Mean  Squared  Error  Versus  the  Number  of  Iterations  for  Three 
Additive  White  Noise  Variances. 

Additionally,  the  learning  time  (convergence  time)  is  not  appreciably  affected  by  the 
noise  contamination  level. 

8.5.B.  Random  Range  Motion  Effects 

The  effects  of  correlated  random  range  motion  on  the  squared  error  are  shown 
in  Figure  8.3.  Correlated  range  errors  are  of  interest  because  they  model  the  errors 
that  might  occur  with  a  range  tracking  control  loop.  The  range  errors  will  be 
correlated  because  the  radar  has  a  finite  servo  bandwidth.  The  simulated  range  errors 
are  modeled  by  a  Gaussian  noise  vector  that  has  been  low  pass  filtered  with  a  second 
order  Butterworth  filter  [5].  In  the  test,  the  learning  rate  is  fixed  at  a  value  of ;/  =  0.1. 
For  this  series  of  tests,  the  variance  of  the  range  errors  is  fixed  at  0. 1  squared  range 
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Range  Motion  Effects 


Figvire  8.3.  Mean  Squared  Error  Versus  the  Number  of  Iterations  for  Various 
Range  Tracker  Servo  Bandwidths  and  Fixed  Range  Tracker  Error 
Variance. 

resolution  cells.  The  plot  shows  the  mean  square  error  for  three  values  of  bandwidth: 

B  =  0.01,  0.0316,  0.1  (5  dB  increments),  and  1.0  (uncorrelated).  It  should  be  noted 
that  the  decorrelation  times  for  each  of  the  bandwidths  are  100,  32,  10,  and  1 
iterations,  respectively.  Additionally,  since  the  gradient  descent  technique  converges 
within  approximately  10  iterations,  the  technique  converges  before  the  range  errors 
decorrelate  for  the  two  lower  bandwidth  simulations  (B  =  0.01  and  B  =  0.0316).  From 
the  curves  in  Figure  8.3,  the  mean  square  error  appears  to  increase  as  20  dB/decade 
with  servo  bandwidth  for  the  two  lower  bandwidth  cases.  This  means  that  the 
accuracy  of  the  range  profile  estimates  is  dependent  upon  the  bandwidth  of  the  random 
motion.  Additionally,  the  20  dB/decade  change  suggests  that  the  error  is  strongly 
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related  to  iteration-to-iteration  variation  of  the  motion  (the  time  derivative  of  the  range 
motion  error). 

Figure  8.4  shows  the  effect  of  varying  the  variance  of  the  range  tracker  error 
while  holding  the  servo  bandwidth  of  the  tracker  constant.  For  this  test,  the 
normalized  bandwidth  is  held  at  a  constant  value  of  B  =  0.0316.  The  variance  of  the 
range  tracking  error  is  set  to  values  of  0.0316,  0.1,  and  0.316  (5  dB  increments). 

From  the  plot,  it  can  be  seen  that  the  mean  square  error  of  the  range  profile  estimate 
increases  as  10  dB/decade  with  increasing  range  tracker  error  variance.  Since  the 
range  profile  estimation  error  increases  20  dB/decade  with  servo  bandwidth  and  10 
dB/decade  with  range  tracker  error  variance,  it  is  reasonable  to  assume  that  the  mean 
square  error  of  the  range  profile  estimate  is  dependent  upon  the  derivative  of  the  range 
tracker  error  signal  and  not  the  range  tracking  error  itself. 

The  effect  of  the  learning  rate  on  the  mean  square  range  profile  error  in  a 
range  motion  environment  is  shown  in  Figure  8.5.  As  is  shown  in  Figure  8.5,  as  long 
as  the  learning  parameter  is  greater  than  the  bandwidth,  the  mean  square  error 
decreases  with  increasing  learning  rate.  In  other  words,  the  learning  rate  should  be  set 
so  that  the  gradient  descent  technique  converges  faster  than  the  range  tracker  errors 
decorrelate.  Additionally,  it  appears  that  the  mean  square  range  profile  error  decreases 
as  20  dB/decade  with  increasing  learning  rate. 

8.5.C.  Random  Angle  Motion  Effects 
In  this  section,  the  effect  of  angular  scintillation  on  the  gradient  descent 
method  is  investigated.  From  Chapter  n,  angular  motion  has  two  effects  on  the  range 
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Range  Motion  Effects 


Figure  8.4.  Mean  Squared  Error  Versus  the  Number  of  Iterations  for  Three  Values 
of  Range  Tracker  Error  Variance  with  a  Fixed  Servo  Bandwidth. 


Range  Motion  Effects 


Figure  8.5.  Mean  Squared  Error  Versus  the  Number  of  Iterations  for  Three  Values 
of  Learning  Rate  and  Fixed  Range  Tracker  Error  Variance. 
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profile:  a  phase  modulation  effect  and  an  envelope  effect.  From  Equation  (2.21),  the 
phase  modulation  varies  significantly  for  a  15  meter  length  target  at  approximately  2 
GHz  if  the  angle  varies  by  more  than  0.14°.  Additionally,  from  Equation  (2.22)  with 
AR  =  0.5  meters,  the  range  envelope  varies  significantly  if  the  angle  varies  by  any 
more  than  1.92°.  As  in  the  previous  section,  angle  motion  is  modeled  by  a  second 
order  Butterworth  low  pass  filtered  Gaussian  random  process.  The  variance  of  the 
process  is  held  constant  to  0.1  squared  degrees  (standard  deviation  of  0.316°).  The 
curves  of  Figure  8.6  represent  averages  of  100  independent  simulations  with  the 
bandwidth  B  =  0.01,  0.0316,  0.1  (5  dB  steps),  and  1.0  (uncorrelated).  The  range 
profiles  are  generated  by  Xpatchl  at  5°  below  waterline  in  elevation  and  from  5°  to 
15°  in  azimuth  sampled  every  0.1°  to  satisfy  Equation  (2.23).  The  range  profiles  for 
azimuth  angles  between  the  calculated  points  are  determined  by  linear  interpolation  of 
the  calculated  points. 

Figure  8.6  shows  that  the  mean  square  error  of  the  range  profile  estimation 
increases  approximately  as  20  dB/decade  with  the  bandwidth  of  the  angle  motion 
when  the  range  errors  do  not  decorrelate  within  a  learning  period  of  approximately  10 
iterations.  This  is  similar  to  the  phenomenon  seen  in  Figure  8.4  for  the  range  tracking 
errors.  Therefore,  as  with  the  range  tracking  errors,  it  appears  that  the  rate  of  change 
in  the  target  incident  angle  is  the  most  important  factor  in  range  profile  estimation. 

Figure  8.7  shows  that  the  mean  square  error  of  the  range  profile  estimation 
increases  as  10  dB/decade  with  the  variance  of  the  angle  motion  (assuming  fixed 
bandwidth).  This  trend,  along  with  those  of  Figure  8.6,  suggests  that  the  rate  of 
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Angle  Motion  Effects 


Figure  8.6.  Mean  Squared  Error  Versus  the  Number  of  Iterations  for  Various 
Angle  Motion  Bandwidths  and  Fixed  Motion  Variance. 


Angle  Motion  Effects 


Figure  8.7.  Mean  Squared  Error  Versus  the  Number  of  Iterations  for  Three  Values 
of  Angle  Motion  Variance  and  Fixed  Motion  Bandwidth. 
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change  of  the  target  incident  angle  is  more  important  that  the  variance  of  the  process. 
Again,  this  is  identical  to  the  error  sources  associated  with  range  tracking  errors. 

Figure  8.8  shows  the  effect  of  changing  the  learning  rate  while  maintaining  a 
fixed  motion  bandwidth  of  B=0.0316  and  an  angular  motion  variance  of  0.1  squared 
degrees.  As  above  with  the  range  errors,  the  higher  the  learning  rate,  the  lower  the 
mean  square  range  profile  error.  This  also  suggests  that  the  source  of  error  is  in  the 
rate  of  change  of  the  range  profile  due  to  angular  motion. 

Angle  Motion  Effects 


Figure  8.8.  Mean  Square  Error  Versus  the  Number  of  Iterations  for  Three  Values  of 
Learning  Rate  and  Fixed  Motion  Variance. 

8.6.  Conclusions 

In  this  chapter,  an  iterative  signal  design  technique  was  derived  and  simulated. 
The  derivation  was  based  on  the  well  known  method  of  gradient  descent.  The 
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computer  simulations  were  used  to  develop  an  understanding  of  the  performance  of  the 
algorithm  in  order  to  determine  a  possible  design  philosophy  for  the  technique. 

From  the  simulations  it  was  seen  that  there  are  basically  two  types  of  errors 
that  are  associated  with  the  iterative  technique:  additive  noise  errors  and  target  profile 
fluctuation  errors.  The  additive  noise  errors  are  representative  of  additive  white  noise, 
jamming,  and  ground  clutter  contamination.  The  target  profile  fluctuation  errors  result 
from  the  fact  that  the  target  range  profile  is  itself  a  random  process.  The  target  is  a 
random  process  since  the  range  to  the  target  and  the  aspect  angle  of  the  target  are 
constantly  changing  by  unknown  amounts. 

The  iterative  algorithm  handles  the  two  types  of  errors  in  different  ways. 
Additive  errors  are  minimized  by  reducing  the  learning  rate.  In  other  words,  the  lower 
learning  rate  causes  the  additive  noise  components  to  be  averaged  out.  Target  range 
profile  fluctuation  errors  are  reduced  by  increasing  the  learning  rate.  This  allows  the 
iterative  algorithm  to  effectively  track  the  range  profile.  If  the  learning  rate  is  set  too 
low,  the  estimated  target  range  profile  converges  to  an  average  range  profile  because 
of  the  random  fluctuations.  Therefore,  a  simple  algorithm  design  rule  would  be  to  set 
the  learning  parameter  in  such  a  way  as  to  mirumize  the  errors  due  to  the  additive 
noise  components  without  increasing  the  errors  due  to  target  range  profile  fluctuation. 
This,  of  course,  will  depend  on  the  level  of  additive  noise  contamination  and  the 
fluctuation  rate  of  the  target.  If  the  target  is  slowly  fluctuating,  then  the  target  can  be 
extracted  from  a  harsh  additive  noise  environment.  If,  on  the  other  hand,  the  target  is 
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fluctuating  rapidly,  only  a  minimal  amount  of  additive  noise  removal  can  be 
accomplished. 

It  should  also  be  pointed  out  that  the  performance  metric  (the  mean  square 
error),  is  calculated  based  on  the  present  value  of  the  range  profile,  while  the 
estimated  range  profile  represents  the  average  range  profile.  Therefore,  the 
performance  metric  is  biased  toward  non-fluctuating  targets.  The  actual  representation 
of  the  target  range  profile  may  be  better  than  the  performance  metric  implies. 

In  this  chapter,  only  one  target  was  used  for  the  simulations  (Target  1  from 
Chapter  VII).  From  other  simulations,  not  shown  in  this  study,  the  algorithm  appears 
to  perform  similarly  on  various  targets. 
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CHAPTER  IX 


SUMMARY  AND  CONCLUSIONS 

The  purpose  of  this  study  was  to  develop  a  signal  processing  scheme  that 
would  combine  the  separate  research  efforts  of  low  probability  of  intercept  signaling, 
low  altitude  low  observable  target  tracking,  and  noncooperative  target  recognition  into 
a  single  research  effort.  The  end  result  of  this  fusion  is  a  new  radar  design  concept 
designated  adaptive  FSK/PSK  signaling.  An  adaptive  FSK/PSK  radar  system  makes 
use  of  matched  FSK/PSK  signals  in  conjunction  with  a  gradient  descent  adaptive 
signal  design  algorithm,  matched  FSK/PSK  signals  are  ultra  wideband  pseudorandom 
multiple  frequency  bi-phase  modulated  signal  sequences  that  have  spectral  content 
specifically  matched  to  the  target  under  track.  This  use  of  target-specific  information 
in  the  design  of  radar  signals  is  the  salient  point  of  this  research  effort.  It  allows  for 
more  efficient  use  of  signal  power  and  results  in  a  performance  improvement  over 
existing  radar  technologies  that  are  matched  to  the  transmit  signal  only.  Since  the 
spectral  content  of  the  target  is  not  known  a  priori,  the  target  range  profile  must  be 
learned  by  the  radar.  This  matching  of  the  signal  to  the  target  is  accomplished 
through  the  use  of  an  iterative  (adaptive)  algorithm  that  is  based  on  gradient  descent. 
Once  the  radar  has  adapted  to  the  target  (learned  the  target),  the  range  profile  of  the 
target  is  known  and  target  identification  can  take  place  based  on  the  learned  range 
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profile.  Therefore,  an  adaptive  FSK/PSK  radar  system  can  be  thought  of  as  a  learning 
system  as  opposed  to  a  measurement  system.  In  other  words,  the  radar  measurement 
is  improved  as  a  result  of  the  learning  process. 

9.1.  Conclusions 

The  pseudorandom  sequence  structure  of  the  matched  FSK/PSK  signal,  along 
with  its  high  time-bandwidth  product,  makes  the  radar  inherently  difficult  to  detect. 
Therefore,  the  signals  can  be  classified  as  low  probability  of  intercept.  Additionally, 
since  the  transmitter  power  is  used  in  a  most  effective  manner  (signal  efficiency  is 
maximized  by  target  matching),  lower  peak  power  signals  can  be  utilized  by  the  radar. 
This  results  in  an  additional  improvement  in  LPI  performance. 

In  an  additive  noise  environment,  matched  FSK7PSK  signaling  provides  a 
signal  processing  gain  over  a  baseline  radar  that  utilizes  a  signal  with  an  impulse-like 
ambiguity  function.  This  performance  improvement  results  from  the  fact  that  the 
signal  power  is  utilized  in  an  efficient  manner  and  the  target  reflection  is  enhanced 
without  enhancing  the  additive  noise.  The  level  of  signal  processing  gain  is  dependent 
upon  the  spectral  complexity  of  the  target.  If  the  target  reflects  aU  frequencies  equally 
well,  little  processing  gain  is  seen.  On  the  other  hand,  if  the  target  reflectivity  is 
concentrated  in  a  few  isolated  frequency  bands,  then  the  processing  gain  can  be  large 
(on  the  order  of  the  square  of  the  number  of  frequencies  used  in  the  signal). 

In  a  ground  clutter  limited  environment,  matched  FSK/PSK  signaling  performs 
better  than  conventional  radar  technologies  that  do  not  take  into  account  any  target 
specific  information.  This  occurs  because  the  matched  FSK/PSK  radar  system  is 
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capable  of  achieving  a  theoretical  target  enhancement  of  between  1  and  while  the 
integrated  sidelobe  level  is  increased  by  a  factor  of  between  1  and  K  where  K  is  the 
number  of  frequencies  used  in  the  construction  of  the  signal.  Therefore,  in  a  clutter 
limited  environment,  adaptive  FSK/PSK  radar  systems  can  achieve  performance 
improvement  factors  of  between  1  and  K.  It  should  be  noted  that  a  matched  FSK/PSK 
radar  system  and  a  conventional  radar  system  using  a  matched  filter  post-processor  are 
able  to  achieve  identical  performance  in  a  clutter  limited  environment. 

Analysis  of  predicted  real  target  signature  data  shows  that  the  use  of  target- 
specific  information  in  the  design  of  the  radar  signal  typically  yields  average  signal  to 
noise  ratio  improvements  of  between  3  dB  and  6  dB  depending  on  the  target.  In  a 
clutter  limited  environment,  typical  improvements  are  between  2  dB  and  5  dB.  These 
improvements  are  closer  to  the  low  end  of  the  theoretically  possible  performance 
improvements.  This  is  result  of  the  uniform  nature  of  optical  region  target  spectra. 

The  gains,  however,  are  significant  and  provide  a  significant  performance  increase.  It 
should  be  pointed  out  that  in  the  resonance  region,  greater  gains  are  expected  since 
resonance  region  spectra  are  typically  impulse-like. 

The  adaptation  algorithm  derived  in  Chapter  VIH  provides  a  method  of  not 
only  designing  the  optimum  signal  to  excite  the  maximum  response  from  the  target, 
but  also  for  measuring  the  best  mean  squared  estimate  of  the  target  range  profile. 
Through  computer  simulations,  it  can  be  seen  that  the  adaptation  algorithm  removes 
additive  noise.  The  level  of  noise  removal  is  dependent  on  the  value  of  the  learning 
parameter.  The  smaller  the  learning  parameter  (slower  the  learning  rate),  the  better  the 
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mean  squared  estimate.  The  simulations  also  show  that  the  random  fluctuations  of  the 
target  range  profile  due  to  random  range  motion  and  random  incidence  angle  motion 
increase  the  mean  squared  error  of  the  estimate.  If  the  learning  rate  is  high  enough 
for  the  adaptive  algorithm  to  learn  the  target  range  profile  before  the  profile 
decorrelates,  then  the  errors  are  small.  If  on  the  other  hand,  the  target  decorrelates 
more  rapidly  than  the  algorithm  can  learn,  range  profile  estimation  errors  are 
increased.  This  results  in  decreased  target  enhancement  due  to  matched  FSK/PSK 
signaling.  It  should  be  pointed  out  that  the  performance  metric  used  in  this  study  (the 
mean  squared  range  profile  error)  is  biased  towards  stationary  (nonfluctuating)  targets. 
It  is  difficult  to  ascertain  the  deleterious  effects  of  range  profile  fluctuations  on  an 
actual  target  identification  algorithm.  The  effects  may  be  insignificant  since  the 
adaptation  algorithm  determines  a  best  mean  squared  fit  to  the  second  order  statistics 
of  the  random  target  fluctuation  process  within  the  learning  period  of  the  algorithm.  It 
is  of  interest  to  note  that  a  linear  averaging  scheme  would  yield  a  zero  result  for  range 
profile  estimation  since  the  phase  angle  of  the  returns  is  random  from  measurement  to 
measurement. 

In  conclusion,  adaptive  FSK/PSK  signaling  radar  is  an  excellent  platform  for 
fulfilling  the  future  mission  of  radar;  low  probability  of  intercept  noncooperative  target 
recognition.  It  utilizes  low  probability  of  intercept  signals  in  a  spectrally  adaptive 
manner.  The  ability  of  the  radar  to  adapt  its  spectrum  makes  it  possible  to  enhance 
targets  thus  improving  system  performance  in  additive  noise  environments  and  ground 
clutter  limited  environments.  Finally,  the  adaptation  scheme  derived  in  this  study  is 
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capable  of  extracting  the  range  profile  of  the  target  under  track.  This  is  important  for 
target  identification  since  it  is  the  range  profile  that  is  the  feature  vector  used  by  most 
target  identification  algorithms. 

9.2.  Future  Work 

The  purpose  of  this  study  was  to  fuse  three  existing  research  efforts  into  a 
single  new  integrated  effort.  Therefore,  sigmficant  work  remains  to  be  done  in  this 
new  research  area.  Three  questions  are  of  particular  interest.  First,  there  are  several 
techniques  for  optimizing  functions  other  than  gradient  descent  [1],  It  would  therefore 
be  of  interest  to  investigate  other  optimization  techniques  within  the  framework  of 
adaptive  FSK/PSK  radar.  Second,  algorithms  exist  for  optimizing  the  spanning  phase 
sequence  of  a  signal  in  order  to  minimize  a  known  clutter  distribution  [2].  It  is  of 
interest  to  determine  if  theses  clutter  minimization  techniques  could  be  applied  to 
adaptive  FSK/PSK.  Third,  it  is  of  interest  to  determine  the  allowable  target  variation 
in  range  and  velocity  versus  processing  time  requirements  to  determine  the  allowable 
performance  envelope  for  the  iterative  radar. 

It  is  also  of  interest  to  point  out  some  other  areas  of  applicability  of  matched 
FSK/PSK  signaling  that  were  not  addressed  specifically  in  this  study.  This  work  can 
be  applied  to  any  area  of  active  sensor  signaling  such  as  underwater  sonar  for  close- 
range  mine  detection,  ultrasound  diagnostics  in  the  field  of  medicine,  fault  detection  in 
the  inspection  of  composite  materials  in  heavy  industry,  and  also  in  the  area  of  laser 
radar.  Additionally,  matched  FSK/PSK  signaling  could  be  utilized  in  the  field  of 
spread  spectrum  communications  through  harsh  fading  channels.  This  application 


would  require  that  the  transmitter  transmit  only  in  spectral  regions  that  are  not 
severely  attenuated  thus  improving  the  signaling  efficiency. 
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CAR_hebb . m 

B.  J.  Skinner 
5  November  1994 
Mississippi  State  University 

This  program  simulates  the  gradient  descent  signal 
design  technique  derived  in  Chapter  VII.  The  range 
profiles  used  in  the  simulation  are  generated  by 
xpatchl,  and  the  matlab  script  read_data_ll .m  must 
be  run  to  read  in  the  xpatchl  frequency  domain  data. 
The  matlab  script  trans_data.m  must  then  be  run  to 
transform  the  data  into  the  time  domain.  The  data 
files  are  aspect  dependent  in  azimuth.  Interpolation 
is  used  to  determine  the  range  profile  for  angle  data 
that  is  not  predicted  by  xpatch.  The  azimuth  angle 
variation  is  modeled  as  a  2nd  order  Butterworth 
low  pass  random  process.  The  range  to  the  target 
is  also  modeled  using  a  2nd  order  Butterworth  low 
pass  random  process  to  simulate  a  range  tracking 
system  with  random  errors. 


This  portion  of  the  script  reads  in  all  needed  data  from 
the  user. 

Mbig=2*K; 

k=menu (' SELECT  LEARNING  RATE' , 'mu=l' , 'mu=. 5' . . 

' mu= . 1 ' , ' mu= .05',' mu= .01'); 

if  k==l 

mu=l; 

elseif  k==2 
mu= . 5 ; 
elseif  k==3 
mu= . 1 ; 
elseif  k==4 
mu=. 05; 
elseif  k==5 
mu= .01; 

end 

k=menu ( ' SELECT  SNR' , ' -6  dB' , ' -3  dB' , ' 0  dB'  , . 

'3  dB' , ' 6  dB' , ' 9  dB' , ' inf inite' ) ; 

var=sqrt (2 . ^ (-k+3) /K) ; 

if  k==7 
var=0 ; 

end 

hl=zeros (K, 1) ; 

k=menu(' SELECT  THE  NUMBER  OF  ITERATIONS' ,' 50' . . 
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'  100' , '  200' , ' 400' , '  800' ) ; 
Nit=50*2" (k-1) ; 


k=menu (' SELECT  TRAINING  EPOCHS' ,' 1' ,' 10' . . 

' 100' , ' 1000' ) ; 

Nave=10  . ''  (k-1)  ; 

clear  m_recl 
clear  m_rec2 
clear  m_rec3 

for  epoch=l:Nave 

%For  each  averaging  realization,  a  new  angle  track  is 
%created. 

angle_bw_s=0 .0316; 
angle_sigma_s=l . 779; 

[b_az,  a_az] =butter (2, angle_bw_s) ; 

az_pntrl=f iltfilt (b_az, a_az, angle_sigma_s* . 

randn (1, 8*Nit) ) +beam_sector; 
ll=5*Nit; 

%For  each  averaging  realization,  a  new  range  track  is 
%created. 

rng_bw= .0316; 
rng_sigma=0 . 0; 

[b_rng,  a_rng] =butter (2, rng_bw) ; 

rng_pntrl=filtfilt (b_rng, a_rng, rng_sigma* randn (1, 8*Nit) ) ; 
12=5*Nit; 


%Re-initialize  all  used  vectors. 

chi= [zeros (1, K/2-1)  1  zeros (1, K/2) ]' ; 
az=az (1, : ) ; 
score=zeros (l,Nit) ; 
shift=zeros (l,Nit) ; 
avging=zeros (size (chi) ) ; 
avg_prof=zeros (size (chi) )  ; 

%Clear  required  memory, 
clear  h2 
clear  reel 
clear  rec2 

clear  corr_score_chi 
clear  corr_score_avg 
ft  clear  rec_pwr 

clear  min_score 

^  %Interpolation  of  the  angle  data. 

%This  should  be  used  only  when  the  angle  does 
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%not  vary  with  the  simulation.  Otherwise,  comment 
%it  out  as  shown  below.  It  should  only  be  done 
%to  speed  up  the  program, 

%  for  m=l:K 

%  h2 (m) =interpl (az, real (11 (m, : ) ) , az_pntrl (11+1)  - 

%  , '  linear' ) + j*interpl (az, imag (11 (m, :))... . 

%  ,  az_pntrl (11+1) ,' linear' ) ; 

%  end 

%  h_prof=h2'; 


for  i=l:Nit 

11=11+1; 

12=12+1; 


%Calculate  angle  interpolation.  Should  be 
%commented  out  when  angle  does  not  vary  within 
%the  simulation.  Most  time  consuming  part  of 
%the  code, 
clear  h2 
for  m=l:K 

h2 (m) =interpl (az, real (11 (m, :))... . 

, az_pntrl (11+1) , ' linear' ) . . . . 

+ j*interpl (az, imag (11 (m, :))... . 

, az_pntrl (11+1) , ' linear' )  ; 

end 

h_prof=h2' ; 

%Calculate  frequency  domain  modulation  for  range  shifting. 
%This  is  done  by  first  calculating  the  range  dependent 
%linear  phase  shift  in  the  frequency  domain.  This 
%Modulation  is  then  multiplied  by  the  FFT  of  the  target 
%range  profile,  and  the  product  is  inverse-FFT  back  into 
%the  range  domain.  This  makes  it  possible  to  shift  the 
%range  profile  a  fraction  of  a  resolution  cell. 

phase2=exp  ( j*2*pi* (fo/B+ (0:  (K-1) ) /K) *rng_pntrl (12+1)  )  ; 
h2=ifft (fft (h_prof ) . *phase2' ) ; 


%Renormalize  the  range  profile. 
h2=h2/sqrt (h2' *h2) ; 

%Add  in  white  Gaussian  noise. 

noise2= (var) /sqrt (2) * . 

(randn (size (h2) ) + j*randn (size (h2) )  )  ; 

%Gradient  descent  design  method. 

chi=chi+2*mu* (h2+noise2) * (h2+noise2) ' *chi; 
chi=chi/sqrt (chi' *chi) ; 
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clc 

%Calculate  performance  Metrics. 

corr_score_chi (i)=max(abs (ifft (fft (abs (chi) 

conj (fft (abs (ref_prof) ))))); 

[corr_score_chi (i)  corr_score_avg (i) ] 
perf=Ti  2* (l-corr_score_chi (i) ) ] 

end 


%Plot  range  track  erros. 
figure (1) 
subplot (2,1,1) 

plot (1 : Nit , rng_pntrl (5*Nit+l : 12) ) 
axis (  [1  Nit  -1  1] ) 
grid  on 
subplot (2, 1,2) 

plot (1 ;Nit, az_pntrl (5*Nit+l : 12) ) 
axis([l  Nit  beam_sector-l  beam_sector+l ] ) 
grid  on 


%Plot  the  learning  curve 
figure (2) 
subplot  (2,1,1) 
plot  (1 :  64,  abs  (chi)  ) 
axis  (  [  1  64  0  1 ] ) 
subplot  (2, 1,2) 

loglog  (1  :Nit,  2*  (l-corr_score_chi)  ,  1  :Nit,  rec_pwr.  ''2) 
axis([l  Nit  .0001  1]) 

if  epoch==l 

m_rec2=2* (l-corr_score_chi)  ; 
m_rec3=rec_pwr .  ^'2; 
else 

m_rec2= [m_rec2;  2* (l-corr_score_chi) ]  ; 
m_rec3=  [m_rec3;  rec_pwr .  2  ]  ; 
figure  (3) 

loglog (1 : Nit, mean (m_rec2) , 1: Nit, mean (m_rec3) ) 
axis([l  Nit  .0001  1]) 
end 

end 
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%  r  e  a  d_cia  t  a_l  i  n .  m 

%  This  program  is  used  to  read  in  an  Xpatchl  data 
%  and  form  it  into  a  usable  Matlab  matrix.  All  four 
%  polarizations  are  retained  (W,  VH,  HV,  and  HH)  . 

%  The  Xpatch  output  data  must  be  pre-processed  by 
%  the  FORTRAN  program  named  field_mat . for 

%Read  in  the  data  file. 
m_file=input (' Enter  the  .lin  filename; 
m_f ile_i= [m_f ile  ' . lin' ] ; 
fid=fopen (m_file_i, '  r'  )  ; 
disp (' Reading  the  datafile') 

[N] =f scanf (f id, ' %d  %d  %d',[3  1]) 

X=fscanf  (fid, '%g  %g  %g  %g  %g  %g  %g  %g  %g  %g  %g' .  [H  inf] 
fclose (fid) 

%Strip  out  the  data  into  separate  vectors. 

disp (' Stripping' ) 

el=X(l,  :)  ; 

az=X (2, : ) ; 

f=X(3,:); 

LL=X(4, : ) +j*X(5, : ) ; 

LR=X(6,  : )  +j*X(7,  :)  ; 

RL=X(8,  :)+j*X(9,  :)  ; 

RR=X(10,  :)+j*X(ll,  :)  ; 

clear  X 

N_freq=N (1) ; 

N_el=N (2) ; 

N_az=N (3) ; 

%Reshape  the  data  vectors  into  matrices 
disp ( ' Reshaping' ) 
el=reshape (el, N_f req, N_el*N_az)  ; 
az=reshape (az, N_f req, N_el*N_az) ; 
f=f (1 :N_f req)  ; 

LL=reshape (LL,N_freq,N_el*N_az)  ; 

LR=reshape (LR, N_f req, N_el*N_az)  ; 

RL=reshape (RL,N_freq,N_el*N_az) ; 

RR=reshape (RR, N_freq, N_el*N_az)  ; 
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%trans_data.m 

%  This  program  converts  the  frequency  domain  data 
%  generated  by  read_data_lin  into  time  domain  range 
%  profile  vectors. 

%Transform  the  data  into  the  time  domain, 
disp (' Transforming' ) 
for  i=l : N_el*N_az 

11  ( : , i) =f ft shift (ifft (LL ( : ,  i)  )  )  ; 

Ir { : , i) =f ft shift (if ft (LR ( : ,  i)  )  )  ; 
rl  ( : , i) =fftshift (ifft (RL ( : ,  i)  )  )  ; 
rr  ( : ,  i)  =fft shift  (ifft  (RR( : ,  i)  )  )  ; 
end 

%Normalize  the  data  to  unit  energy. 

pll=ones  (N_freq,  1)  *  (std(ll)  .  (-1)  )  /sqrt  (N_freq)  ; 

plr=ones  (N_freq,  1)  *  (std(lr)  .  (-1)  )  /sqrt  (N_freq)  ; 

prl=ones (N_freq, 1) * (std(rl) (-1) ) /sqrt (N_freq) ; 

prr=ones  (N_freq,  1)  *  (std(rr)  . (-1)  )  /sqrt  (N_freq)  ; 

ll=pll.*ll; 

lr=plr . *lr; 

rl=prl . *rl; 

rr=prr . *rr; 


o\o  o\o  o'? 


params .m 

This  matlab  script  contains  the  initial  constants 
required  by  the  CAR_hebb.iu  simulation  script. 
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B=300e6; 
N=1024; 
K=64; 
fo=2 . 0e9; 
c=3e8; 

T=K/B; 

Td=N*T; 


beam  sector=10; 


131 


%This  matlab  script  reads  in  a  data  file  containing 
%radar  cross  section  data  and  calculates  the  target 
%enhancement ,  signal  to  clutter  improvement,  and 
%LPI  performance  loss  due  to  Matched  FSK/PSK  signaling. 
%This  program  is  used  to  calculate  the  data  contained 
%in  Table  7.1. 

m_file=input ('Enter  the  .lin  filename:  ','s'); 
m_f ile_i= [m_file  '.lin']; 
fid=fopen (m_file_i,  '  r' )  ; 
disp (' Reading  the  datafile') 

[N] =f scanf (f id, ' %d  %d  %d',[3  1]) 

X=fscanf (fid, ' %g  %g  %g  %g  %g  %g  %g  %g  %g  %g  %g',[ll  inf]); 
fclose (fid) 

disp (' Stripping' ) 
el=X (1, : ) ; 
az=X (2, ; ) ; 
f=X(3, : ) ; 

LL=X(4,  : )  +j*X(5,  :)  ; 

LR=X(6,  :)+j*X(7,  :)  ; 

RL=X(8,  : )  +j*X(9,  : )  ; 

RR=X (10, : ) + j*X (11, : ) ; 

clear  X 

N_f req=N (1) ; 

N_el=N (2) ; 

N_a2=N (3) ; 

disp ( ' Reshaping'  ) 

el=reshape (el, N_f req, N_el*N_az) ; 
az=reshape (az,N_freq,N_el*N_az) ; 
f=f (1 :N_freq)  ; 

LL=reshape (LL,N_freq,N_el*N_az) ; 

LR=reshape (LR, N_f req, N_el*N_az) ; 

RL=reshape (RL,N_freq,N_el*N_az) ; 

RR=reshape (RR, N_f req, N_el*N_az) ; 

num=sum (abs (LL) .^2) ; 
denl=max (abs (ifft (LL) ) )  ; 
den2=sum (abs (LL)  )  ; 

Eb= (num. /deni . /den2)  .^2; 

Em=N_f req*num. /den2; 
sc=num.  /  (deni .  ''2)  /N_freq; 

L=l/N_freq. /max (abs (LL) . / (ones (N_freq, 1) *den2) ) ; 

perf= [10*logl0 (mean (Eb) )  10*logl0 (max (Eb) ) ; . 

10*logl0 (mean (Em) )  10*logl0 (max (Em) )  ; . 


10*logl0 (mean (sc) )  10*logl0 (max (sc)  ) 
10*logl0 (mean (L) )  10*logl0 (min (L)  )  ] 
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C  FIELD_MAT.F 

C  THIS  FORTRAN  PROGRAM  READS  IN  AN  XPATCHl  DATA  FILE 

*  C  AND  CONVERTS  THE  DATA  INTO  A  FORM  THAT  IS  EASILY 
C  PROCESSED  BY  THE  MATLAB  SCRIPT  READ_DATA_LIN . M. 

C  ADDITIONALLY,  THE  VERTICAL  POLARIZED  DATA  IS  CONVERTED 
'*■  C  INTO  CIRCULAR  POLARIZED  DATA  IF  NEEDED. 

COMPLEX  W,  VH,  HV,  HH 
COMPLEX  RR,  RL,  LR,  LL 
COMPLEX  J 

REAL  F,  AZ,  EL 

INTEGER  N_FREQ,  N_AZ,  N_EL,  N_POINTS 
CHARACTER* 50  JUNK 

CHARACTER  F_FILE*13,  M_FILE*13,  L_FILE*13 
CHARACTER* 5  FILE 

Q-k  *  *  *  *  -k  -k  *******  -k  *******************  -k  *****  ir  -k  *****  -k  *********  * 

C  GET  THE  REQUIRED  DESCRIPTION  FROM  THE  USER 

Q*  *********************************  -k  **********************  * 

PRINT  *, 'ENTER  THE  TARGET  FILENAME . FIELD' 

READ  (*,*)  FILE 

* 

PRINT  *, 'ENTER  THE  NUMBER  OF  FREQUENCIES' 

READ  (*,*)  N_FREQ 

% 

PRINT  *, 'ENTER  THE  NUMBER  OF  AZIMUTH  ANGLES' 

READ  ( * , * )  N_AZ 

PRINT  *,' ENTER  THE  NUMBER  OF  ELEVATION  ANGLES' 

READ  ( * ,  * )  N_EL 

N_POINTS=N_FREQ* (N_AZ+1) * (N_EL+1) 

F_FILE=FILE  //  '.field' 

M_FILE=FILE  //  '.cps' 

L_FILE=FILE  //  '.lin' 

print  *,F_FILE,  M_FILE,  L_FILE 

PRINT  *,  " 

PRINT 

*''THE  FILE  TO  BE  PROCESSED  IS  ',F_FILE 

PRINT  *,'THE  OUTPUT  FILE  IS 
PRINT  *,'THE  NUMBER  OF  FREQUENCIES  IS 
PRINT  *,'THE  NUMBER  OF  AZIMUTH  ANGLES  IS 
*  PRINT  *,'THE  NUMBER  OF  ELEVATION  ANGLES  IS 

PRINT  *,'THE  NUMBER  OF  DATA  POINTS  IS 
PRINT 


=='  PRINT 

'  ,M_FILE 
'  , N_FREQ 
'  ,N_AZ  +  1 
'  ,N_EL+1 
' ,N  POINTS 
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Q********************************************************* 

C  OPEN  FILES  FOR  PROCESSING 

*  Q*****************************************************-**** 

OPEN(UNIT=l,  FILE=F_FILE) 

OPEN(UNIT=2,  FILE=M_FILE) 

OPEN(UNIT=3,  FILE=L_FILE) 

Q********************************************************* 

C  READ  THE  DATA,  PROCESS  IT  AND  WRITE  IT  OUT  * 

Q********************************************************* 

J=CMP LX  (0.000, .1.000) 

WRITE (2,*)  N_FREQ,  N_AZ+1,  N_EL+1 
WRITE (3,*)  N_FREQ,  N_AZ+1,  N_EL+1 
DO  1  ICOUNT=l,5 
1  READd,*)  JUNK 

ICOUNT=0 

10  READ  (1,  *,  END=100)  EL,  AZ,  EL,  AZ,  F,  W,  VH,  HV,  HH 
ICOUNT=ICOUNT+l 

LL=.5*  (  (HH-J*VH)  -J*  (HV-J*W)  ) 

LR= .  5*  (  (HH-J*VH)  +  J*  (HV-J*W)  ) 

•  RL=.  5*  (  (HH+J*VH) -J*  {HV+J*W)  ) 

RR=.5*  (  (HH+J*VH)+J*  (HV+J*W)  ) 

>  WRITE (2,*)  EL,  AZ,  F,  REAL (LL) ,  AIMAG(LL), 

1  REAL (LR) ,  AIMAG (LR) , 

■  2  REAL(RL),  AIMAG (RL)  , 

3  REAL(RR),  AIMAG (RR) 

WRITE  (3,*)  EL,  AZ,  F,  REAL  (W)  ,  AIMAG  (W)  , 

1  REAL (VH) ,  AIMAG (VH) , 

2  REAL (HV) ,  AIMAG (HV) , 

3  REAL(HH),  AIMAG (HH) 


GOTO  10 

100  PRINT  *,'ALL  DONE' 

PRINT  *,'YAKNOW .  THERE  WERE  ',ICOUNT,'  DATAPOINTS' 


* 


END 


